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ABSTRACT 


Energetic  chemical  species  have  been  investigated  by  both  theoretical 
and  experimental  methods.  The  nature  of  the  chemical  bonding  in  H3  has 
been  studied  by  means  of  a  natural  spin  orbital  analysis  of  a  previously 
calculated  wavefunction.  H3+  is  show  to  resemble  closely  its  united  atom 
analog,  Li,  and  the  Hartree-fock  energy  of  H3+  is  estimated  to  be  -1.301 
hartrees.  Also,  electron  density  plots  are  presented  in  order  to  clarify 
the  physical  picture  of  the  bonding  in  H3  ,  and  suggestions  are  made  as  to 
how  further  calculations  might  proceed  most  efficiently.  In  the  geminal 
theory  of  chemical  bonding  a  new  method  has  been  developed  for  solving  the 
coupled  integro-differential  equations  for  determination  of  the  optical 
gsminals  for  a  system.  Geminal  theory  has  also  been  generalized  to  include 
odd  electron  systems.  A  microwave  spectrometer  employing  Zeeman  modulation 
and  a  flow  cell  has  been  built  for  studying  transient  paramagnetic  species. 
Charge  storage  in  yirradiated  polymethylmethacrylate  (PMMA)  at  room 
temperature  has  been  studied  by  measuring  thermoelectric  currents.  By 
using  pure  PIMA  as  well  as  PMMA  doped  with  hole  and  electron  traps  it  has 
been  shown  that  added  electron  traps  but  not  hole  traps  lead  to  charge 
stabilization.  The  thermoelectric  currents  appear  to  be  due  to  hole 
mobility.  The  trapped  charge  in  doped  PMMA  decays  with  an  activation 
energy  of  0.9  ev.  Polyvinyl  chloride  stabilizes  even  larger  amounts  of 
charge  than  does  doped  PMMA. 
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STRUCTURE  AND  REACTIVITY  OF  ENERGETIC  CHEMICAL  SPECIES 


SCOPE  OF  THE  PROGRAM 

The  physics  and  chemistry  of  highly  energetic  chemical  species  such  as 
ions,  electronically  and  vibrat ionally  excited  molecules,  and  free  radicals 
is  still  in  its  infancy.  The  greatest  paucity  of  information  concerns 
excited  state  species;  their  structure,  reactivity,  and  metastability 
represent  a  newly  developing,  experimentally  difficult,  but  potentially 
rewarding  field  of  study.  Theoretical  techniques  are  just  now  arriving  at 
tne  point  at  which  they  can  provide  accurate  and  detailed  information  on 
energetic  species  of  chemical  interest  (e.g.  H3+).  Thus,  theory  is 
reaching  the  point  where  it  can  supplement  and  guide  the  experimental  in¬ 
vestigation  of  energetic  species.  The  characteristics  of  energetic  species 
of  all  types  in  the  solid  state  is  also  an  area  of  great  importance  and 
limited  information. 

The  areas  of  study  are  vectored  toward  an  enhanced  understanding  of 
the  physical  and  chemical  characteristics  of  highly  energetic  chemical 
species.  A  balanced  experimental  and  theoretical  study  of  both  structure 
and  reactivity  of  carefully  selected  species  and  systems  utilizing  a 
variety  of  techniques  character i zes  the  general  approach  used.  In  this 
annual  report  progress  is  reported  in  the  following  areas:  (a)  nature 
of  bonding  in  H3+,  (b)  geminal  theory,  (c)  construction  of  a  microwave 
spectrometer  for  studying  transient  energetic  species,  and  (d)  radiation- 
induced  charge  storage  in  organic  polymers. 


I.  THEORETICAL  INVESTIGATIONS 


1.  INTRODUCTION 

The  present  theoretical  aspect  of  the  program  is  directed  toward  the 
theoretical  elucidation  of  the  structure,  reactivity,  and  general  nature 
of  highly  energetic  species.  In  order  to  provide  accurate  and  detailed 
information  concerning  compounds  of  potential  interest  as  a  rocket  pro¬ 
pellant  as  soon  as  possible,  a  multi-faceted  program  has  been  begun  in 
which  both  specific  compounds  and  general  techniques  are  being  studied 
simultaneously. 

The  vork  designed  to  elucidate  the  nature  of  particular  species  was 
centered  on  the  H3+  molecule,  both  because  it  holds  considerable  interest 
as  a  reactive  species,  and  because  it  can  be  used  as  a  prototype  for  the 
description  of  the  bonding  in  electron  deficient  systems. 

The  other  area  in  which  considerable  progress  has  been  made  is  in 
the  development  and  implementation  of  a  geminal  theory  of  chemical  bonding. 
This  theory  can  be  applied  generally  to  give  accurate  results  on  systems 
of  chemical  interest,  without  losing  all  hope  of  interpreting  the  resulting 
wavefunction.  Since  geminal  theory  is  still  in  its  infancy,  much  funda¬ 
mental  analysis  was  necessary  before  computational  algorithms  could  be 
developed,  and  much  of  this  analysis  was  completed  during  this  year.  In 
particular,  significant  advances  have  been  made  in  two  separate  aspects  of 
pair  theory.  First  of  all,  a  new  method  has  been  developed  for  solving 
the  coupled  integro-differential  equations  that  arise  in  the  determination 
of  the  optimal  geminals  for  a  system.  This  method 


appears  to  hold 
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considerable  promise  toward  significantly  reducing  the  computational  time 
that  is  necessary  to  determine  the  best  set  of  geminals,  without  loss  of 
accuracy.  Secondly,  it  has  been  feasible  to  show  that  it  is  not  only 
possible,  but  convenient,  to  include  an  odd  electron  into  the  general 
theory  of  geminals,  so  that  the  concept  of  geminals  need  not  be  restricted 
to  systems  with  an  even  number  of  electrons. 

As  well  as  the  aforementioned  developments  in  the  area  of  analysis, 
much  of  the  work  necessary  to  implement  ab  initio  calculations  using 
geminal  theory  has  been  carried  out.  Most  of  the  computer  programs  to 
evaluate  the  one-  and  two-center  integrals  that  arise  in  these  calculations 
have  now  been  written  and  tested  on  the  GE-625  computer  at  the  University 
of  Kansas.  This  work  involved  a  careful  choice  of  the  manner  in  which 
these  programs  were  implemented.  The  result  has  been  that  the  integral 
programs  are  presently  very  efficient  and  often  run  3  to  4  times  faster  than 
other  existing  programs  without  loss  of  accuracy. 
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2.  ANALYSIS  Of  H3+ 
a.  Purpose 

Although  the  small,  reactive  species,  H3+,  has  been  known  to  exist  for 
a  considerable  length  of  time,  the  lack  of  experimental  characterization  of 
H3+  beyond  confirmation  that  it  exis'.s  has  required  that  most  knowledge  of 
this  molecule  must  come  from  theoretical  considerations.  To  investigate 
the  nature  of  Hj*,  the  following  analysis  of  an  accurate  wavefunction  was 
carried  out. 

The  use  of  a  wavefunction  which  is  derived  from  an  arbitrarily  chosen 

basis  for  interpretation  of  the  nature  of  the  bonding  in  two-electron 

systems  often  encounters  severe  difficulties.  In  a  series  of  papers, 

it  has  been  emphasized  that  if  one  is  to  attempt  to  extract  the  real  'nature' 

of  the  bond,  then  one  must  make  use  of  some  invariant  properties  which  are 

not  solely  a  function  of  the  arbitrary  basis  set  chosen.  There  are  various 

such  invariant  descriptions  of  a  molecular  system.  The  most  commonly  used 

7 

one  is  the  Hartree-Fock  approximation,  or  a  close  approximation  thereto. 

This  particular  invariant  description  is  an  important  one  because  energy  is 

used  as  a  criterion  of  excellence  and  because  it  can  be  shown  that  most 

one-electron  properties  are  given  correctly  to  second  order  by  the  exact 

Hartree-Fock  function.  Calculations  now  exist  for  an  extensive  series  of 
8 

molecules,  and  informative  and  valuable  correlations  have  been  made  using 
these  results. 

Unfortunately,  the  Hartree-Fock  energy  is  far  from  the  exact  energy, 
usually  farther  than  the  entire  binding  energy  of  the  molecule.  Conse¬ 
quently,  sublet ies  in  the  molecular  structure  are  imperfectly  revealed  by 
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studies  of  Hartree-Fock  functions  alone.  Also,  higher  invariant  approxima¬ 
tions  using  minimum  energy  as  a  criterion  of  excellence  rapidly  become 
extraordinarily  difficult  to  organize.^ 

As  an  alternative  invariant  description,  the  eigenfunctions  and  eigen¬ 
values  of  the  first  order  density  matrix  associated  with  a  given  trial 
function  have  been  used.  These  eigenvalues  are  the  occupation  numbers  of 
the  natural  orbitals,  and  even  for  quite  approximate  functions,  have  been 
found  to  be  almost  basis  independent10  in  the  systems  examined  thus  far. 
Consequently,  conclusions  that  can  be  drawn  from  the  occupation  numbers  of 
natural  orbital  expansions  appear  to  have  a  range  of  validity  that  is  both 
esthetically  pleasing  and  heuristically  useful. 

The  results  of  a  natural  orbital  analysis  of  the  H?+  wavef unction 
reported  by  Christoffersen11  are  presented  here.  Since  the  energy  of  this 
wavefunction  is  close  to  that  expected  from  the  exact  wavef unct ion,  it  is 
likely  that  it  contains  representative  components  of  all  important  basis 
types.  It  appears  that  this  analysis  is  not  significantly  different  from 
one  that  would  be  obtained  from  the  exact  wavefunction,  or  from  any  other 
computed  approximate  wavefunction  that  has  as  low  an  energy  expectation 
value  as  the  wavefunction  used. 

b.  The  Wavefunction 

We  limit  our  consideration  to  the  equilibrium  position  of  the 
ground  A)  state  of  .  "We  have  previously  discussed  in  another  connection 
natural  orbitals  and  occupation  numbers  associated  with  some  very  crude 
functions  for  the  linear  conformation^  but  as  Christoffersen  has 
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demonstrated,  the  equilateral  conformation  is  by  far  the  most  stable.  The 
wavef unction  analyzed  here  is  of  the  form 

■b  =  nT  (1)  c  n  (2),  (i) 

in  which  n  is  a  column  vector  of  normalized  (but  not  necessarily  orthogonal) 
~  and  qT  |s  the  transpose  of  q. 

symmetry  adapted  orbitals,/  In  Table  l-l  we  list  the  precise  definitions 

of  the  individual  q  orbitals  in  terms  of  the  original  basis  set  of  Slater 

type  orbitals  on  each  nucleus.  Here  a.  A,  a,  an  are  centered  on  the  first 

vertex;  b,  B,  2,  or  jin  on  the  second  vertex,  and  c,  C,  y,  yn  on  the  third 

vertex  of  the  triangle.  Lower  case  Roman  letters  refer  to  Is  orbitals: 

l  _ 

(a,  b,  c)  =  (C ,  /n)2  e“  ir ,  (2) 

capital  Roman  letters  to  2s  orbitals: 

(A,  B,  C)  -  (C*  /3")5  re"^2r,  13) 

lower  case  Greek  letters  to  2p  orbitals 

5 

C  -  -  c 

(a,  *1,  y)  =  (-~)*  r  cos  8  e  3  .  (4) 

The  p  orbitals,  a,  y,  are  directed  from  each  nucleus  towards  the 
midpoint  of  the  molecule.  The  p  orbitals,  an,  j}n,  yn  are  directed  perpen¬ 
dicular  to  the  plane  of  the  molecule.  Chr i stoff ersen* s  calculations 
showed  that  p  orbitals  in  the  plane  perpendicular  to  a,  (3,  y  respectively 
were  relatively  unimportant,  and  these  were  omitted  from  the  calculation 
analyzed  here.  The  internuclear  distance  was  taken  at  R  =  1.6575  bohrs 
and  the  parameters  chosen  to  be  C,  =  1.20,  =  1.175,  Cj  =  1*80,  This 

function  gave  a  computed  energy  of  -1.33264  hartrees. 
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Table  l-l.  Normalized  symmetry  orbitals.  See  text 
for  definition  of  original  basis. 
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In  Table  l-ll  we  exhibit  the  coefficient  matrix  £.  It  should  be  noted 
that  these  coefficients  are  different  from  those  reported  by  Christ offer sen 
because  we  have  now  normalized  the  individual  molecular  orbitals  composed  of 
Slater-type  basis  orbitals.  Since  the  original  basis  is  chosen  symmetry- 
adapted,  the  complete  £  matrix  has  many  zeros  from  symmetry.  Only  the 
non-zero  diagonal  blocks  are  tabulated  in  Table  l-ll.  Because  of  the 
doubly  degenerate  irreducible  representation  in  the  pertinent  group  D.^, 
there  is  a  duplication  of  blocks  as  indicated  in  the  table.  The  3-3,  6-6, 
and  9-9  configurations  were  omitted  in  the  Christoffersen  wavefunct ion. 

It  has  already  been  indicated  that  the  original  basis  is  non-orthogonal. 
In  Table  l-ll I  we  tabulate  the  non-zero  overlap  integrals  between  the 
normalized  basis  functions  f).  One  of  the  striking  features  of  these 

numbers  is  the  high  values  of  many  of  the  overlaps,  especially  that  of  A12 

connecting  the  totally  symmetric  Is  and  2s  functions.  These  values  make 
questionable  any  attempts  to  make  use  of  the  coefficients  of  Table  l-ll  in 
obtaining  a  handle  on  the  real  nature  of  the  wavefunct ion.  It  is 

precisely  for  this  reason  that  one  strongly  prefers  to  discuss  the  wave- 

function  in  terms  of  its  natural  expansion. 
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Table  1-11.  The  coefficient  matrix  £  of  equation  1. 


lisas* 

Value 

Element 

Value 

*11 

1.385331 

*44  "  *77 

-0. 123776 

Cl*  ■  C*1 

-O.285686 

*45  “  *54  "  *78  "  *  87 

0. 016589 

c13  “  C31 

0.092328 

C46  “  *84  “  *79  “  *97 

-0.022103 

cat 

0.008291 

*55  *  *88 

-0.008559 

**»  *  csa 

0.001684 

*58  ”  *85  “  *89  “  C96 

0.006762 

cs3 

0 

*88  “  *99 

0 

*10,10 

-0.037462 

*11,11  “  *12, IS 

-0.006999 

All  other  elementa  are  identically  zero  by  symnetry. 
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Table  l-lll.  Overlap  matrix  A  over  the  original  baaia. 


Aw  -  -  0.910431 

A13  -  A31  -  0.487879 

Asa  ■  Ais  •  0.276857 


4*5  -  ^54  -  -  0.712956 

4*e  ■  4b4  "  ^79  ■  A^  •  -0.777956 
^e  "  A*  “  4bs  *  ^ae  ■  -0.742063 


All  diagonal  elements  of  A*  1,  All  other  elements  of  A  that 
are  not  explicitly  given  are  zero. 
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c.  The  Natural  Expansion 


Since  our  baai*  i«  non -orthogonal,  it  ia  convenient  to  uae  the  non -orthogonal 
for  obtaining  the  natural  orbitala  aa  outlined  by  Barnett,  Linderberg 


and  Shull.12  The  proceaa  ia  straightforward.  Since  we  have  a  two  electron  problem. 
It  is  eaaier  to  diagonalise  C  A  rather  than  J_  A .  The  reaulta  are,  of  courae, 
equivalent  aa  abown  by  L^wdln  and  Shull.15  Me  obtain  natural  orbitala  of  ai.  s’, 
a£,  and  e"  aymmetry.  Since  there  waa  only  oue  molecular  orbital  uaed  of  a£  type, 
and  one  degenerate  pair  of  e"  type,  the  natural  orbital  occupation  nuabera  of 
these  are  slsguy  given  by  the  squares  of  the  correapondlng  coefficienta  c10,i0. 

*«•  cia.ia  of  Tabl*  11 •  0,1  tb*  other  h*nd*  5x5  “trlx  Product,» 

■uet  be  diagonal  lead  to  obtain  the  natural  orbitala  of  ax  and  e'  symmetry.  The 
reaulta  are  Hated  in  Table  IV.  The  natural  expansion  of  the  wavef unction  la 

then  given  by 

(5) 

tad  the  correapondlng  expansion  of  the  first  order  density  matrix  by 


with  omission,  of  course,  of  the  usual  spin  factors  correapondlng  to  a  pure 
singlet  function. 

y,  consider  first  the  occupation  numbers.  To  place  these  in  proper  context 

+ 

V9  coolidtr  those  of  the  united  atocn,  Li  > 

and  those  of  the  singlet  coupled  separated  atoms  H  +  H  +  H  ,  each  an  infinite 
distance  from  the  others.  For  the  latter,  the  correct  uavefunctlon  for  the 

ground  state  la  given  by  ^ 

V  \  a(l)b(2)  +  b(l)a(2)  +  a(l)c{2)  +  c(l)a(2)  +  b(l)c(2)  +  c(l)b(2)j  . 
omitting  the  usual  singlet  spin  function.  The  occupation  numbers  for  this  case 
are  2/5  for  the  natural  orbital. (a  +  b  +  c),  and  ±  each  for  a  degenerate  pair 
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which  might  be,  for  example,  (a  -  b)  and  (b  -  c)  or  any  independent 
linear  combinations  of  these.  The  s^nmetry  of  the  first  approximate 
natural  orbital  in  this  case  is  and  the  degenerate  pair  is  e' ,  just 
as  is  found  for  the  actual  H3+  function  as  analyzed  in  Table  l-IV.  On  the 
other  hand  we  note  that  the  H3*  equilibrium  state  function  has  a  much 
higher  a,  occupation  number,  and  correspondingly  lower  values  for  e' . 

The  situation  is  more  complex  for  the  united  atom.  Natural  orbitals 
in  this  case  are  of  the  usual  types,  s,  p,  d,  etc.  One  needs  to  know  the 
decomposition  of  these  degenerate  spherical  harmonics  into  terms  appropriate 
for  the  lower  symmetry  Djf,.  Obviously  s  becomes  a', .  The  three  p  levels 

N  j  |  f| 

split  into  a2  ♦  e' ;  the  five  d  levels  split  into  a,  *  e  ♦  e  ;  and  so  on. 

One  can  make  the  very  crude  approximation  that,  even  though  the 
spherical  degeneracies  have  been  split,  the  occupation  numbers  remain 
equally  divided  among  the  components.  This  is  expected  to  be  a  good 

approximation  only  at  the  united  atom,  but  it  allows  a  zeroth  order  com- 

10 

par  ison  of  occupation  numbers.  Shull  earlier  has  used  a  similar  division 
of  correlation  energies  in  comparing  He  and  H2  with  moderate  success.  Thus 
we  assign  one-third  of  the  total  p  occupation  in  Li  to  a2  and  two-thirds 
to  the  two  components  of  e' .  We  use  occupation  number  data  secured  by 
Macias1^  from  the  35  configuration  function  for  Li*  by  ieiss.  ^  The 
first  six  terms  with  total  occupation  number  0.999971  may  thus  be  divided 
into  0.998121a',,  0.000594a2,  0.001223e',  and  0. 000033a". 

The  total  aj  occupation  number  (n,  *  n3  ♦  n12)  of  0.985680  for  Hj* 
thus  falls  nicely  into  the  sequence  with  the  Li+  a,  total  of  0.996121  and 
the  separated  atom  total  of  0.666667.  The  united  atom  character  of 
is  clearly  evidenced.  Correspondingly,  the  e'  occupation  total  has  climbed 
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Table  I  — I V .  Natural  Orbitals  for  H«. 


Syaaatry 

C1 

ni 

K 

■ 

*1 

0. 992404 

0.984865 

*. 

-0.079810 

0. 006370 

X3 

a’ 

-0.079810 

0. 006370 

Xa 

it 

*2 

-0.037462 

0.001403 

1 

*1 

-0.028547 

0.000815 

Xa 

a" 

-0.006999 

0.000049 

a" 

-0.006999 

0.000049 

X. 

a1 

-r.  006107 

0.000037 

X. 

a1 

-0.006107 

0.000037 

Xjlo 

a' 

+0.001592 

0. 000002s 

An 

a' 

+0.001592 

0.000002s 

%xa 

1 

•l 

+0.000592 

O.OOOOOO4 

1.000000 

1.191445  m  -  0.275794  7  2  +  0.093129  Is 
-1.215053  q4  +  0.081936  ns  -  0.220262  ne 
-1.215053  nr  +  0.081936  ne  -  0.220262  n» 

n  10 

0.978720  ni  -  1.566924  n2  +  0.573740  n3 

Hu 
n  12 

1.02543674  -  1.246616  n5  +  0.355811  7e 
1.025436  nr  -  1. 246616  n  e  ^  0-355811  7  e 
0.546234  1  4  +  0.960382  n  S  +  1-708073  7  a 
0.546234  n  7  ♦  0. 960382  n  e  +  1.708073  7  0 
2.576951  7  1  -  2.216009-12  *  1.154724  ns 


u 

from  the  0.001223  of  the  united  atom  to  0.012819  in  H3+  towards  O.333333 
for  the  separated  atoms.  On  the  other  hand,  the  occupation  number  of 
0.000594  in  the  united  atom  increases  to  O.OGI4O3  in  H3+  and  subsequently 

disappears  in  the  separated  atoms.  The  same  kind  of  maximum  is  observed 
for  the  e"  occupation. 

The  occupation  number  data  are  therefore  not  unexpectedly  different 
from  the  united  atom  case,  tending  toward  the  separated  atoms.  The  main 
changes  are  in  the  direction  of  transferring  occupation  from  s  type  to 
p  type  terms  associated  with  the  transfer  of  in-out  type  of  correlation  to 
a  more  formally  angular-type  correlation  which  is  more  advantageous  with 
separated  nuclear  centers.  This  data  is  summarized  in  Table  l-V. 

d.  Truncated  Expansions 

Several  calculations  have  been  reported  recently  in  which  the  Hartree- 
Fock  energy  of  equilateral  ^  has  been  estimated.  Using  a  Gaussian  basis 
set  of  s-  and  p-type  orbitals,  Lester  and  Krauss16  report  a  Hartree-Fock 
energy  of  -1.27778  hartrees  at  R  =  1.68  bohrs.  Preuss17  used  a  larger 
Gaussian  basis  set  of  s-type  orbitals  and  obtained  a  Hartree-Fock  energy  of 
-1.2988  hartrees  at  R  =  1.62  bohrs.  Also,  in  a  one-center  calculation  using 
s,  d,  f,  g  and  h  orbitals,  Joshi18  reported  a  Hartree-Fock  energy  of 
-1.2863  hartrees  at  R  =  1.2523  bohrs.  Although  a  Hartree-Fock  calculation 
was  not  explicitly  carried  out  by  the  present  authors,  it  is  possible  to 
obtain  a  reliable  estimate  of  the  Hartree-Fock  energy  from  the  analysis  of 
the  wavefunction  reported  here. 

It  has  been  shown  that  the  energy  calculated  from  a  truncated  natural 
expansion  of  one  term  (renormalized)  is  very  close  to  the  Hartree-Fock 
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onefgy,  the  error  being  of  fourth  order.  ^ *  Since  the  wavefunction  used 

in  this  analysis  is  quite  flexible,  it  is  expected  that  the  energy  of  the 

first  natural  orbital  will  be  quite  close  to  the  exact  Hartree-Fock  energy. 

The  energy  obtained  using  the  first  natural  orbital  was  -1.2992  hartrees 

^  ~  1 *6575  bohrs,  and  possesses  an  overlap  with  the  final  wavefunction 

of  0.9024038.  Since  the  optimum  SCF  energy  has  been  observed  to  occur  at  a 

slightly  shorter  distance  than  the  observed  equilibrium  distance3  in  H2, 

it  might  be  expected  that  the  value  of  -1.2992  hartrees  might  be  lowered 

slightly  by  performing  a  calculation  at  the  equilibrium  distance  observed 

17  18 

in  other  SCF  calculations  *  (R  =  1.62  bohrs).  However,  it  has  also 

been  observed  that  the  potential  surface  is  quite  flat  with  respect  to 

1? 

variation  of  R.  On  this  basis  it  was  anticipated  that  the  energy  im¬ 
provement  obtained  by  such  a  calculation  would  be  negligible,  and  it  was  not 
carried  out. 

The  proximity  of  the  present  value  of  -1.2992  hartrees  to  the  Hartree- 

Fock  energy  can  be  further  elucidated  in  the  following  manner.  For  a 

similar  two-electron  system,  HeH+,  S.  Peyerimhoff  found  that  a  basis  set 

consisting  of  only  six  basis  orbitals  (o1su  ,  o2su  ,  <J2S'  ,  o2p,,  ,  ols  ,  o2d  ) 

He  He  He  He  H  rH 

with  optimized  orbital  exponents  gives  a  Hartree-Fock  energy  within  approxi¬ 
mately  0.0008  hartrees  of  the  true  Hartree-Fock  energy.  Since  the  basis  set 
of  Is,  2s,  2po,  2pn  orbitals  employed  for  the  calculations11  on  H3+  is  very 
similar  to  that  used  by  Peyerimhoff,  including  optimization  of  orbital 
exponents,  it  seems  likely  that  the  error  in  the  Hartree-Fock  energy  of 
H3+  is  probably  only  slightly  larger  than  the  error  in  the  six  basis  orbital 
calculation  of  Peyerimhoff.  The  error  due  to  the  additional  nucleus  in 
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H3  is  not  easily  estimable,  but  if  we  allow  0.001  hartrees  as  the  error 
estimate  of  this  complication,  we  may  combine  this  with  the  error  estimate 
of  0.0008  hartrees  associated  with  a  similar  basis  set,  and  estimate  the 
Hartree-Fock  energy  of  H3*  to  be  approximately  -1.301  hartrees.  A  surwnary 
of  recent  Hartree-Fock  calculations  is  given  in  Table  |— VI . 

each  of  the  two  electron  species  that  is  closely  related  to  H3+  has  a 
correlation  energy  which  is  approximately  constant,  with  H2,  HeH+,  and  Li* 
having  the  values  0.041  hartrees,22  0.042  hartrees21  and  O.O43  hartrees,23 
respectively.  Since  the  orbital  shapes  differ  greatly  in  these  three 
species,  it  appears  that  the  shape  of  an  orbital  does  not  play  a 
significant  role  in  determining  the  correlation  energy  of  two-electron 
systems.  It  has  also  been  shown24  that  the  correlation  energy  for  two-electron 
atoms  remains  essentially  constant  at  O.O4O-O.O45  hartrees,  regardless  of 
the  nuclear  charge.  Since  the  size  of  a  given  orbital  is  determined  largely 
by  the  nuclear  charge,  the  constancy  of  the  correlation  energy  for  iso- 
electronic  two  electron  atoms  implies  that  the  size  of  the  orbital  is  also 
an  unimportant  factor  in  determining  the  correlation  energy.  Since  the 
orbitals  of  H3+  differ  from  the  orbitals  in  the  isoelectronic  analogs 
moat  1 y  in  size  and  shape,  it  is  likely  that  the  correlation  energy  for  H  ♦ 
will  also  lie  within  the  range  C.O4O-O.O45  hartrees.  If  we  take  0.042 
hartrees  as  the  estimated  correlation  energy  for  H/,  then  we  obtain  a 
value  of  -1.343  hartrees  as  the  estimated  true  total  energy  of  H3*. 

It  is  interesting  to  note  that  this  estimate  of  the  total  energy  is 
considerably  above  the  value  of  -1.357  hartrees  reported  by  H.  Conroy.25 
Since  the  method  used  by  Conroy  is  not  constrained  by  the  variation  theorem, 
it  is  possible  that  the  extrapolated  results  may  lie  sither  above  or  below 
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Table  l-V. 

Summary  of  Occupation 

Number  Data 

_ 

e' 

a  " 
a2 

e" 

Li+ 

0.998121 

0.001223 

0.000594 

0.000033 

H,+ 

0.985680 

0.012819 

0.001403 

0.000098 

2H  +  H+ 

0.666667 

0.333333 

0.000000 

0.000000 

Table  I -VI,  Summary  of  Hartree-Fock  Calculations  on  Equilateral  H3+ 


Reference 

Method 

R^(bohrs) 

ehf 

(hartreea) 

W.  A.  Lester  and 

M.  Krauss* 

Gaussian  orbitals 

1.68 

-1.27778 

B.  D.  Joshi^ 

One-center  expansion. 
Slater-type  orbitals 

1.623 

-1.2863 

H.  PreussC 

Gaussian  orbitals 

1.62 

-1.2988 

M.  E.  Schwartz** 

Gaussian  orbitals 

1.6505 

-1 .2999 

Present  Results 

First  natural  orbital 

1.6575 

-1.2992 

Estimated  "true" 
Hartree-Fock  energy 

- 

- 

-1.301 

a  Reference  16. 
k  Reference  18, 
c  Reference  17. 


**  Private  communication. 
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the  true  energy.  In  a  similar  comparison  of  calculations  on  HeH  ,  Michels 

pointed  out  that  the  extrapolated  energy  obtained  by  Conroy  appeared  to  be 

too  low,  since  the  extrapolated  value  fell  below  that  of  a  generalized  64- 

27 

term  James-Cool idge  type  wavefunction  used  by  Wolniewicz.  Although  the 
wavefunction  employed  by  Chr istoffersen  for  H3+  (12  configurations)  is  not 
as  flexible  as  the  one  used  by  Wolniewicz  for  HeH  ,  it  is  probably  as 
flexible  as  the  11  configuration  wavefunction  for  HeH+  used  by  Anex,^  which 
diffeis  by  only  0.0037  hartrees  in  total  energy  from  the  Wolniewicz 
wavefunction.  Furthermore,  the  use  of  a  Slater-type  basis  set  for  Hj  by 
McClean,  Weiss  and  Yoshimine  similar  to  that  used  For  H3+  resulted  in  an 
error  of  0.007  hartrees  in  the  total  energy.  Consequently,  it  appears  that 
the  error  of  -0.01  hartrees  to  be  added  to  the  total  energy  of  -1.33264 
hartrees  from  the  H3+  calculation11  as  implied  above  is  a  reasonable 
estimate,  and  that  the  value  of  -1.357  hartrees  reported  by  Conroy  may  be 
too  low. 

Since  the  natural  orbital  expansion  of  rank  two  can  always  be 
expressed  in  "u,v''  form  which  has  the  physical  intuitional  simplicity  of 
the  Hei tier -London  approach,  it  is  also  of  interest  to  find  the  energy  of 
the  natural  expansion  of  rank  two  for  the  present  wavefunction.  In  this 
case,  the  "two  term"  expansion  is  actually  composed  of  three  terms, 
necessitated  by  the  degeneracy  of  X2  andX3*  The  energy  obtained  for  this 
truncated  (and  renormalized)  function  was  -1.326  hartrees  at  R  =  1.6575  bohrs. 
The  overlap  of  this  function  with  the  complete  wavefunction  is  0.9988089  at 
the  equilibrium  distance. 
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e.  A  Symmetric  Representation 

Visual  interpretation  of  three-fold  systems  is  very  much  hindered 
by  the  necessary  presence  of  doubly  degenerate  representations,  e'  and  e". 
Although  these  two  representations  as  a  whole  have  three-fold  density,  they 
do  not  provide  an  easy  way  to  discuss  three  equivalent  relatively  localized 
parts  which  may  be  associated  with  (in  this  case)  H-H  "bonds".  For 
example,  the  real  representat ions,  n4  and  n?  in  Table  I -I  do  not  treat  the 
three  atoms  equivalently  in  each  orbital.  One  can  avoid  this  by  going  to 
an  imaginary  representation  with,  for  example, 


e 


3(1  — S ) 


(a  ♦  wb  +  w*c) 


(7) 


and 


= 


1 


3(1-3)  J 


(a  +  w*b  +  wc) 


(8) 


each  of  which  clearly  has  3-fold  symmetry.  But  then  there  is  still  no 
convenient  localized  formulation,  for  the  orbitalsX,»  e{  and  el  are  not 
three-fold  equivalent. 

We  suggest  that  an  attractive  alternative  formulation  can  be  obtained 
by  choosing  one  localized  representative  of  the  e'  pair,  and  generating 
two  equivalent  (but  non-orthogonal)  representatives  by  successive 
application  of  the  Cj  operation.  Thus  consider 


where  Xa>  Xa  form  the  degenerate  pair.  Suppose  y8  to  be  aa 
localized  as  possible  in  one  "bond".  Then  we  chooae  as  bases  for  the 
e'  representation: 

Xr  =  Xa 

XB  =  ca  Xa  "  (-  /T  x3  -  Xs  >/2  (10) 

Xt  =  ^3  Xs  =  (♦  ifT  Xs  -  Xa  )/2. 

We  note  that  these  orbitals  are  normalized  but  are  not 
orthogonal.  Their  mutual  overlaps  are  precisely  -1/2. 

In  this  representation  we  can  now  write 

♦  =  Xi  d)  Xl  <2)  *  |  X*  Cxr  d)  Xr  (2) 

+  X,  d)  X8  <2)  +  Xt  d)  Xt  (2)  ]  , 

and  hence 

♦  =  3  [Xl  Xl  d)  Xi  (2)  -  2  xa  xr  d)  Xr  (2)] 

+  3  CX  J  Xi  d>  Xi  (2)  -  2  X»  X„  d)  X,  U)]  (11) 

+  3  U?  Xi  (1)  Xi  (2)  -  2X,  xt  d)  Xt  (2)]. 

in  which  f  has  now  been  divided  into  three  separate  localized  "sum 

of  squares"  form.  These  may  individually  be  expressed  in  "u,  v" 


form  so  that 
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fa  {[ur  (l)vr(2)  +vr(l)ur(2)]  +  Cu8(l)vi(2)  +  v8(i)ug(2)3 

+  [ut(l)vt(2)  +vt(l)ut(2)]}  /3, 

where 

ur  =  X.i  Xi  +  'ft  ^2  Xr 
vf=  Xi  "  ^  a  Xr> 


(12) 


(13) 


and  similarly  for  the  other  u'  e  and  v's.  At  this  level  of  approximation, 

Hj  can  be  written  as  the  sum  of  three  equivalent  (but  geometrically 
displaced)  two-electron  "bonds"  of  Heitler-London  type. 

It  is  crucial  to  choose  Xs  to  be  as  "bond-like"  as  possible,  i,  e.  , 
most  localised  in  bond  rather  than  on  an  atom.  This  choice  is 
intuitively  clearly  (b  -  c)  /  ^2(1  -  S)\  which  has  no  contribution  from 
atom  a  at  all.  Truncating  the  Christoffersen  wavefunction  to  the  first 
three  natural  orbitals,  and  making  this  phase  choice  for  ya  w®  have 


finally 


uf  a  0.996792  Xi  +  0.  399764  xr 


(14) 


vf  =  0.996792  xi  -  0.  399764  xr- 

Inserting  the  forms  for  Xi  an<*  Xr  ~  Xa  from  Table  IV,  collecting 
coefficients,  we  arrive  at  the  following  representatives  for  u  and  v 
in  the  original  basis: 

ur  a -0.  0715 15*  a  +  0.  464518*  b  +  1.  000551- c 

-0.  0424304*  A  -  0.  0975565-  B  -  0.  152683*  C  (15) 
-0.0269  052*0  +  0.  0442349*  +  0.  115375*  y, 


and 


(16) 
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v  =  1.000551 *a  ♦  0.464518* b  -  0.071515‘c 
-0.1 52683* A  -  0. 0975565* B  -  0. 0424304*0 


+0.115375*a  +  0.0442349*?  -  0.0269052*y, 
where  ur  and  vr  have  not  been  normalized.  Further  visualization  of  localized 
representatives  will  be  given  in  the  following  section. 


f.  Discussion 


We  have  already  seen  that  one  of  the  benefits  of  a  natural  spin 
orbital  analysis  is  an  estimate  of  the  Hartree-fock  energy  of  H3+.  Several 
other  dividends  can  be  gained  via  this  analysis,  not  the  least  of  which  is 
an  indication  of  the  importance  of  various  terms  in  the  configuration 
interaction  expansion. 

One  interesting  point  to  note  concerning  the  occupation  numbers  and 
their  associated  symmetry  types  given  in  Table  I  —  I V  is  that,  with  the 
exception  of  the  aj  representation,  a  natural  orbital  of  each  symmetry  type 
appears  once  before  any  repetition  occurs.  This  suggests  that,  for 
further  configuration  interaction  or  similar  studies  of  H,  ,  and  perhaps 
for  other  molecules  with  similar  bonding  characteristics,  the  most  rapid 
convergence  of  a  many  configuration  wavef unction  will  be  obtained  by  first 
including  molecular  orbitals  of  all  of  the  various  symmetry  types,  rather 
than  adding  further  terms  of  the  same  symmetry.  In  cases  where  only 
limited  size  basis  sets  can  be  employed,  this  principle  should  be  of 
particular  importance.  In  the  case  of  H3+,  the  molecular  orbitals  that 
were  used  represent  all  of  the  possible  symmetry  types  except  a,.  The 
way  in  which  a  molecular  orbital  of  this  type  can  be  formed  is  by  using 
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nucleus  are  chosen  with  the  x,  z  plane  as  the  plane  of  the  molecule,  and 
the  z-axis  pointing  toward  the  center  of  the  molecule.  These  orbitals  give 

•  n  m 

rise  to  two  possible  molecular  orbitals,  belonging  to  the  a,  and  e 
irreducible  representations  of  D^.  Table  I  —  I V  suggests  that  the  use  of 
an  a"  molecular  orbital  in  an  acceptable  configuration  would  result  in  a 
natural  orbital  with  an  occupation  number  greater  than  0.000037,  and 
probably  between  .000049  and  0.000015. 

One  of  the  intriguing  aspects  of  H3+  left  unanswered  in  the  previous 

work^  was  that  of  the  bonding  arrangement  in  In  order  to  elucidate 

the  bonding  in  Hj*,  an  electron  density  plot  was  made  along  a  line  from  one 

nucleus  through  the  center  of  the  molecule,  as  well  as  a  contour  map  of  X2 , 

and  these  plots  are  given  as  Figures  1-1  and  1-2.  As  can  be  seen,  the 

electronic  charge  at  the  center  of  the  molecule  is  greater  than  that  at  the 

midpoint  of  a  line  between  nuclei.  Specifically,  the  value  of  X  at  the 

center  of  the  molecule  is  0.1157,  and  the  value  of  X2  at  the  midpoint 

between  two  nuclei  is  0.1036.  Thus,  if  we  subtract  off  the  large  electron 

density  at  each  nucleus,  there  would  remain  a  build-up  of  density  at  the 

center  of  the  molecule.  This  observation  affords  further  evidence  to 

justify  the  earlier  speculation  that  one  appropriate  description  of  the 

bonding  arrangement  in  Hj+  would  be  a  set  of  three  "bonds"  pointing  from 

the  nuclei  to  the  center  of  the  molecule,  rather  than  by  bonds  along  the 

11 

lines  between  the  nuclei. 

Alternatively,  the  localized  representatives  of  (u,  v)  form  that  were 
introduced  in  the  previous  section  can  be  used  to  describe  the  bonding  if 


Figure  1-3, 


Electron  density  contours  for  the  ground  state  H«  as 
calculated  from  the  first  natural  orbital.  Proceeding 
outward  from  any  given  nucleus,  the  contours  correspond 
to  the  following  values  of  v»:  0.  160,  0.  120,  0.  ll'i,  0.  110 
0.100,  and  0.080.  ' 
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three  equivalent  bonds  are  desired.  However,  these  "bonds"  are 
considerably  "bent",  which  can  be  seen  in  the  following  manner. 

From  equation  (11),  we  write  one  of  the  three  equivalent  bonds* 

as 

*r  (1.2)  =  N  Cur(l)  vr(2)  -r  vr  (1)  ur(2)]  ,  (17) 

We  can  obtain  a  visual  representation  of  by  forming  the  square  of 
♦r  and  integrating  over  the  coordinates  of  electron  two  to  give 

$r  (1)  =  0.284716  [u“  (1)  +  v*  (1)  ]+  0.4H635‘u.(l)vr(l)^  (18) 

which  gives  a  representation  of  the  electron  density  of  electron  one, 

averaged  over  all  positions  of  electron  two.  Contour  plots  of  constant 
a 

$r  (1)  are  given  in  Figure  1-3,  which  clearly  show  the  "bent"  character 

of  the  bond  representation. 

♦ 

Since  H3  can  be  considered  to  be  the  triatomic  analog  of  die 

homonuclear  diatomic  molecule  (H3) ,  it  is  interesting  to  see  if  the 

concepts  of  fractional  natural  alternant  character  and  fractional  natural 

ionic  character  that  have  been  developed  for  homonuclear  diatomics* 

can  be  carried  over  to  the  polyatomic  system.  For  the  diatomic 

a 

systems,  the  fractional  natural  alternant  character  (\^)  and  the  fraction- 

a 

al  natural  ionic  character  (\  +  )  are  defined  by 


3  Oi)*  +  (n8 ) ^  ]  / 

*+3  [{»!)*  -  (n8)*  ]/  VT, 


(19) 


flfM  l-J.  Contour  plots  of  constant  values  of  4^(1).  (See  equation 
(18)  in  the  text. )  The  contours  correspond  to  the  following 
values  of  4$(1):  0. 165,  0. 125,  0. 115  and  0. 10T„ 
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where  n,  and  n8  are  the  occupation  numbers  (renormalized)  of  the  two 
configuration  natural  expansion  of  the  wavefunction. 

In  the  case  of  H3  ,  we  have  already  seen  that  it  is  possible  to 
express  the  molecule  in  terms  of  three  equivalent  "bonds".  However, 
in  order  to  obtain  the  amount  of  alternant  and  ionic  character  for  these, 
"bonds",  it  is  necessary  to  transform  the  expression  describing  these 
"bond."  into  natural  orbital  form.  This  can  easily  be  done  by  construct¬ 
ing  two  new  orthonormal  orbitals,  and  5.,  from  u  and  v  by 


S+  =  (u+v)/N* 
?_  a  (u-v)  /N-  , 

where 


(20) 


Nt  a  Juadv  +  Jv^dV  t  2  JuvdV.  (21) 

The  expression  for  the  "bonds"  in  equation  (17)  can  then  be  written  in 


natural  orbital  form  as 


or 


*  (1,2)  .  N[(EL) 


?♦  -<£■) 


*.  3. 


(22) 


*r(1.2)  =  0.  98731!*  §  +  -0.158800.  5.  . 


(23) 


Using  this  data  in  equation  (19) ,  we  see  that  the  fractional  natural 
alternant  character  and  fractional  natural  ionic  character  of  each  of 
the  three  "bonds"  in  Hs  is  given  by 


a 


0.  6568 


a  0.  3432. 

For  comparison  purposes,  the  same  calculation3  at  the 
equilibrium  di.Unce  (R  -  1. 4  bohrs,  for  H,  result,  in  £  .  0. 394 


(24) 
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and  X2  =  0.606,  and  for  R  =  1.3  bohrs  in  He2+,  the  calculation  results**  in 
X*  =  0.3583  and  X^  =  0.6417.  Thus  we  may  conclude  that  the  degree  of  natural 
ionic  character  of  each  of  the  "bonds"  in  H3+  and  the  bonds  in  H2  and  He*+ 
is  very  similar,  with  H3+  falling  just  below  He2++  and  HU. 

The  observation  concerning  the  nature  of  the  electron  distribution  in 
Hj  is  in  disagreement  with  the  results  of  Preuss.  In  one  set  of  calcula¬ 
tions  using  17  Is-type  Gaussian  orbitals,  he  found  that  the  Hartree-Fock 
energy  of  H3  +  as  calculated  with  two  Gaussians  placed  in  the  center  of  the 
molecule  was  almost  the  same  as  the  energy  obtained  when  the  two  functions 
were  absent  from  the  center.  This  resulted  in  the  conclusion  that  there  was 
only  a  small  electron  density  at  the  center  of  the  molecule.  Unfortunately, 
conclusions  drawn  from  an  examination  of  coefficients  of  terms  involving  a 
non-orthogonal  basis  set  are  vulnerable  to  criticism. 

The  basis  functions  of  Preuss  are  far  from  being  orthonormal,  with 
the  overlap  of  the  various  functions  on  the  nuclei  with  the  two  orbitals 
in  the  center  varying  from  0.21584  to  0.90869.  Consequently,  we  would  like 
to  suggest  that  an  alternative  conclusion  is  possible  from  the  same  data. 

A  calculation  of  the  electron  density  in  the  center  of  the  molecule 
qsing  the  Preuss  wavefunction  gives  \|>2  =  0.118,  compared  to  X2  =  0.116  from 
this  work.  Similarly,  a  point  midway  between  two  nuclei  gives  i|>2  =  0.102, 
compared  to  X2  =  0.104  from  this  work.  Therefore,  although  it  is  not  ap¬ 
parent  from  an  examination  of  the  coefficients  in  the  non-orthogonal 
wavefunction,  the  Preuss  wavefunction  provides  essentially  the  same 
description  of  the  electron  density  in  H3+  as  our  X, .  Furthermore,  the 
correct  form  for  the  electron  density  had  already  been  established  before 
the  addition  of  the  Gaussians  in  the  center  of  the  molecule,  for  a  calculation 
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of  *a,  omitting  <pl6  and  <p1?,  also  results  in  a  larger  density  in  the  center 
of  the  molecule,  compared  to  the  dansity  at  the  midpoint  between  the  nuclei. 
Thus,  the  addition  of  Gaussians  in  the  center  had  little  effect,  but  only 
corrected  small  defects  in  the  electron  distribution  provided  by  the 
Gaussians  already  present.  This  should  not  be  interpreted  as  a  lack  of 
electron  density  at  the  center  of  the  molecule. 

The  above  discussion  might  be  construed  to  imply  that  a  one-center 

calculation  on  Hj  would  be  rather  effect ive*  This  approach  has  met  with 
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only  very  limited  success  *  because,  although  there  is  a  build  up  of 
electron  density  in  the  center  relative  to  the  midpoint  of  a  line  between 
nuclei,  the  electron  density  at  either  of  these  points  is  still  far  lower 
than  at  any  of  the  nuclei.  Therefore,  the  usual  problem  associated  with 
one  center  calculation  is  also  present  in  H}+,  i.e.,  how  does  one  build 
up  charge  density  at  off-center  nuclei?  Thus,  it  is  not  expected  that 
significantly  better  results  will  be  obtained  for  H3+  using  a  one-center 
approach  than  for  any  other  molecular  system  with  off-center  nuclei. 

As  this  work  indicates,  a  great  deal  of  insight  has  been  gained  into 
the  nature  of  H}  ,  but  the  intriguing  question  of  its  lack  of  experimental 
characterization  still  remains.  For  example,  its  extremely  high  binding 
energy  suggests  that  it  ought  to  be  possible  to  stabilize  this  cation  with 
an  appropriate  anion.  Of  course,  one  of  the  crucial  questions  in  this  matter 
is  the  actual  magnitude  of  the  binding  energy  of  H}\  Since  the  work 
described  here  involves  a  variational  calculation,  the  only  information  that 
can  be  gleaned  from  that  calculation  is  that  the  true  binding  energy  is 
larger  tnan  that,  i.e.,  the  calculation  has  given  an  upper  bound  to  the 
true  energy.  Thus,  it  is  of  considerable  interest  to  obtain  a  lower  bound 
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to  tho  true  energ'  of  H3*  to  "bracket"  the  true  energy.  Then  much  more 
definitive  statements  regarding  the  ability  of  H3+  to  combine  with  suitable 
anions  can  be  made. 

The  computation  of  a  lower  bound  usually  is  formulated  so  that  matrix 
elements  over  the  square  of  the  Hamiltonian  for  the  system  are  needed.  In 
the  pasty  the  difficulties  in  evaluating  the  integrals  that  arise  in 
computations  of  this  type  have  made  it  impossible  to  perform  such  computa¬ 
tions  except  for  atoms  and  the  Ha  molecule.  No  method  has  been  developed 
for  polyatomic  molecules.  Howevery  we  have  recently  been  able  to  show  thaty 
if  Gaussian  basis  orbitals  are  employedy  the  integrals  can  be  calculated 
in  closed  form  or  with  only  a  single  numerical  integration.  Fortunatelyy 

a  very  accurate  variational  wavefunction  for  H}+  using  Gaussian  basis 
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orbitals  has  been  given  by  Schwartz.  We  have  succeeded  in  duplicating 
his  original  calculation  for  the  upper  bound  to  Hj+,  and  are  in  the  process 
of  developing  the  necessary  computer  programs  for  the  use  of  the  Schwartz 
wavefunction  to  calculate  a  lower  bound  to  the  ground  state  energy  of  H3*. 

g.  Conclusions 

(1)  The  nature  of  the  electronic  correlation  in  the  ground  state 
of  the  Hj+  molecule  has  been  shown  to  bear  a  striking 
resemblance  to  -  united  atom  analog  (Li  )y  by  means  of  an 
analysis  of  the  occupation  numbers  arising  from  a  natural 
orbital  analysis. 

(2)  The  main  change  involved  in  the  transittor  from  the  united 
atom  to  Hj+  is  that  the  in-out  correlation  in  the  atomic 
case  is  transferred  to  a  more  formally  angular-type 
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correlation*  accompanied  by  a  corresponding  change  in  the 
first  natural  orbital  from  spherical  symmetry  to  the  symmetry 
of  the  Aj  irreducible  representation  of  the  °3h  group. 

(3)  Using  the  calculated  value  of  -1.2992  hartrees  at  H  -  1.6575 
bohrs  for  the  approximate  Hartree-Fock  energy  of  H5  as 
found  from  the  first  natural  orbital,  the  Hartree-Fock  energy 
has  been  estimated  to  be  -1.301  hartrees. 

(4)  Electron  density  plots  of  the  square  of  the  first  natural 
orbital  have  confirmed  the  earlier  speculation  that  the 
appropriate  way  to  describe  the  bonding  arrangement  in  Hj 

is  either  by  means  of  bonds  pointing  toward  the  center  of  the 
molecule  from  each  nucleus,  or  by  means  of  bent  bonds  between 
each  pair  of  nuclei. 

(5)  If  one  assumes  that  the  approximately  constant  value  of  0.042 
hartrees  for  the  correlation  energy  of  other  isoelectronic 
systems  can  be  taken  over  to  ,  then  the  total  energy  of 
H3+  can  be  estimated  to  be  -1.343  hartrees. 
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3.  (SEMINAL  THEORY 

The  work  in  the  area  of  geminal  theory  has  been  focused  towarc  developing 
an  adequate  technique  for  handling  many-electron,  many-center  problems  using 
ab  initio  principles,  but  at  the  same  time  attempting  to  retain  some  degree 
of  intuitive  simplicity  in  the  resulting  wavefunction.  Any  theory  of 
chemical  binding  should  account  for  the  highly  useful  and  successful  experi¬ 
mental  notion  that  the  basic  building  blocks  of  normal  chemical  bonds  are 
pairs  of  electrons,  and  not  individual  electrons  acting  independently.  This 
is  done  in  geminal  theory  by  assuming  that,  for  an  even  number  of  electrons, 
the  total  wavefunction  can  be  written  as  a  properly  ant i symmetrized  product 
of  pair  functions  (geminals); 

*(1,2,...,2N)  =  AG,(1,2)*Ga(3,4)  .  *GN(2N-1 ,2N) , 

where  each  geminal  Gj  is  a  normalized,  antisymmetric  function  of  the  space 
and  spin  coordinates  of  two  electrons.  Furthermore,  each  geminal  can  be  as 
cooplicated  as  the  user  desires,  and  so  can  allow  for  a  very  careful 
description  of  the  interaction  of  two  electrons  within  the  same  geminal. 

It  is  known  e:  perimentally  that,  in  many  cases,  the  interaction  between  two 
electrons  within  the  same  geminal  is  much  stronger  than  the  interaction 
between  electrons  in  different  geminals.  This  is  usually  integrated  into  pair 
theory  via  the  strong  orthogonality  constraint,  i.e., 

Jo’o.UG^I.j)*,  =  6^ 


For  computational  purposes,  each  geminal  is  usually  written  as  a 
product  of  a  space  and  a  spin  part, 

Gu(1,2)  =  Au  (1,2)®  (1,2), 

and  the  space  portion  is  usually  expressed  in  "natural  orbital  form": 

•V 

A  d,2)  =  ^CuiXi(DXi(2), 

i=1 

where  .  is  the  occupation  number  of  the  ith  natural  orbital  in  the  pth 
geminal,  and  is  to  be  determined  so  that  the  energy  is  minimized.  Of 
course,  the  determination  of  the  appropriate  natural  orbitals  Xj  is  also 
done  so  that  the  energy  is  minimized.  Within  these  approximations,  the 

total  energy  of  the  space  portion  of  a  system  consisting  of  2N  electrons 
is  given  by 


an 
geminals  I 


E  = 


^  I  2^Cpi<|1ii hlu i >  *  £  Si iC*j  ia j |iy i] 


all 

geminals 

E 

n  <v 


N  ,i>L 

a  v 


^  .  C2  .  (, 4  [p  .p  .  v  .v  .  -2 fa  •  v  . I  v  -u . 

vjl  J  jJ  f  i  j|  jwi| 


>*J 


where  (ij  ~Xaj>  h  is  a  one-electron  operator,  and 


=  JJ\*(DX?(l)~X^(2)X?)(2)dV,dV2. 


■hen  this  expression  is  minimized  with  respect  to  the  natural  orbital 
occupation  coefficients,  subject  to  the  normalization  condition,  one  must 
•olve  the  following  coupled  eigenvalue  equation  to  determine  C  .  : 

T|K 


i 


•<  =  1,2 . . 

1  =  1,2 . ,  all  geminals. 


"*>*ra  ^ki  *S  4,1  e^Tective  Hamiltonian  for  the  geminal. 

These  equations  are  coupled,  i.e.,  we  need  to  know  all  the  other  C’s  to 
•olve  for  C^.  These  equations  are  usually  easily  solved  via  an  iterative 
technique,  which  begins  by  making  the  initial  approximation  that 


C 


=  i. 


(j>0, 


and  solving  the  above  eqjation  for  £  ,  and  consequently,  a  new  set  of  C  .. 

1  I 

This  procedure  is  then  repeated  until  the  results  are  self-consistent,  i.e. 
the  set  of  Cs  does  not  change  from  one  iteration  to  the  next.  This 
technique  is  not  new,  and  has  proven  to  be  effective  in  practice. 

However,  the  energy  expression  must  also  be  minimized  with  respect  to 
the  natural  orbitals  themselves,  subject  to  the  constraint  that  the  natural 
orbitals  remain  orthonormal,  i.e., 


jX-V  1)Xvj(DdVi  = 

■hen  this  minimization  is  carried  out,  we  obtain  the  following  set  of 

coupled  integro-differential  equations  for  the  determination  of  the  natural 
orbital  Xup 


36 


C2. 

M-  • 


h  + 


£Sn(£jlVj 

jrl 


all 

gemjjnals 

v(*u) 


where 


geminals  Nv 

-  L  L 


V=1  J=1 


9 


JVJ  =  Jx:j(2)XVj(2)  dv,  , 

•<vj‘1>X  (1)  *  { X,j(')X*,j(2)7LX  (2)  , 


and  the  X^.^j  are  Lagrangian  multipliers  which  arise  from  the  orthogonality 
constraints.  Due  primarily  to  the  fact  that  "off-diagonal"  Ls^rar^ian 
multipliers  (X^j^j,  withpi^vj)  appear  in  the  equations,  there  has  not 
to  date  been  a  successful  algorithm  developed  for  the  solution  of  these 
equations.  However,  it  is  shown  below  that  it  is  possible  to  develop  a 
workable  algorithm  for  the  solution  of  these  equations  by  first  transforming 
the  equations  to  a  simpler,  but  entirely  equivalent  form. 

We  first  let 


Mop^> 


V  * 

♦  r  c  .c Vk  .  ♦ 

4-.  HI  UJ  Uj 

J=1 
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geminals 
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and  rewrite  the  equations  determining  the  natural  orbitals  in  matrix  form, 


Mop*  = 


where  X  is  a  column  vector,  given  by 


where  Q  is  the  number  of  geminals,  and 


We  now  note  that  ^  is  a  real,  symmetric  matrix^  This  is  apparent  from  the 
fact  that  anc*  both  arise  from  the  same  orthonormal i ty 

constraint, 


/’CiV'  *  J«*.j 

and  consequently  j ,v j  must  he  identical  to  With  this  information, 


.  •  dV  =  b  *6 .  ■  , 

Pl  'IV  |J  ’ 
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we  note  that  any  real,  symmetric  matrix  can  be  brought  to  diagonal  form  via 
a  unitary  transformation.  Denoting  this  transformation  by  U,  we  can  rewrite 
the  equations  to  be  solved  as: 


or 


=  (uxu+)ux, 


or  M  X* 
op 


u  i 


where  X  =  UxU*  =  diagonal  matrix,  X^.  =£uj,  and  x'  =  UX.  The 
diagonal izat ion  of  tha  original  \  matrix  provides  a  considerable  simplifi¬ 
cation,  because  it  removes  the  off-diagonal  Lagrang’an  multipliers  from 
the  problem,  and  decouples  the  equations  determining  the  natural  orbitals. 

The  solution  must  still  be  done  in  an  iterative  manner,  since  L  (ui) 

op 

depends  upon  a  knowledge  of  all  of  the  other  natural  orbitals.  However, 
the  iterative  nature  of  the  solution  is  entirely  similar  to  that  which 
arose  in  the  determination  of  the  natural  orbital  occupation  coefficients, 

and  it  is  anticipated  that  the  transformation  introduced  here  will 
allow  the  determination  of  the  natural  orbitals  to  be  carried  out  just  as 
easily  as  the  determination  of  the  occupation  coefficients.  The  implementa¬ 
tion  of  the  computer  programs  representing  this  approach  is  in  progress, 
and  a  quantitative  estimate  of  the  time  that  will  be  saved  by  this  analysis 
is  being  sought 

Although  considerable  progress  has  been  mada  toward  development  of  a 
suitable  algorithm  for  solving  the  equations  which  determine  the  geminals, 
gaminal  theory  itself  cannot  be  considered  to  be  a  general  theory  unless 
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the  restriction  to  systems  with  an  even  number  of  electrons  can  be  removed. 
As  shown  below,  it  is  indeed  possible  to  generalize  geminal  theory  to 
include  an  odd  electron. 

Consider  a  system  of  2N+1  electrons,  in  which  2N  of  the  electrons  are 
spin  paired.  The  trial  function  is  written  as  a  properly  antisymmetr ized 
product  of  N  geminals  and  one  orbital: 

*(1,2,  .  .  .  ,  2N,2N+1)  =  A |g,(1,2)G2(3,4)  .  .  .  G^UN-1 ,2N)«(2N*1 ) 

In  a  manner  similar  to  that  in  v^iich  the  individual  geminals  are  treated, 
the  space  portion  of  the  added  orbital  can  be  expanded  in  terms  of  natural 
orbitals,  i.e., 

*(N*D  =  Z  C(N+l)iX(N*1)i  ’ 
i=1 

with 

/^N+l)iK(iVl)jdV  6ij  * 

The  strong orthoganal ity  constraint  is  generalized  to  include 

f  A*(1,2)«(1)dV,  =  0,  |i  =  1,2, . N. 

J  |i 


Wi th  this  form  of  the  wavefunct ion,  the  total  energy  of  the  system  of 
2N+1  electrons  can  be  shown  to  be 
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By  comparison  with  the  previous  energy  expression  for  a  system  in  a  singlet 
state  with  an  even  number  of  electrons*  it  is  seen  that  the  additional 
terms  result  from  the  intergeminal  interactions  caused  by  the  presence  of 
the  additional  orbital. 

As  in  the  singlet  case,  the  energy  expression  must  also  be  minimized 
with  respect  to  both  the  best  natural  orbital  expansion  coefficients  and  the 
natural  orbitals  themselves.  When  this  procedure  is  carried  out  for  the 
doublet  case . equat i ons  analogous  to  those  of  the  singlet  case,  but  slightly 
more  complicated  due  to  the  presence  of  the  additional  intergeminal  terms, 
are  obtained.  However,  the  techniques  that  have  been  developed  for  the 
singlet  case  are  directly  applicable. 
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Thu*  the  analysis  which  is  necessary  to  generalize  geminal  theory  for 
us*  on  systems  of  possible  interest  as  propellants  is  completed,  as  well  as 
a  new  algorithm  for  the  computation  of  the  geminals  which  describe  these 
systems*  Work  is  underway  to  implement  this  analysis  with  the  appropriate 
computer  programs. 


4.  REFERENCES 


1.  H.  Shull,  J.A.C.S. ,  02,  1287  (I960). 

2.  H.  Shull,  J.  Phys.  Chem.,  66,  2320  ( 1962) . 

3.  S.  Hagstrom  and  H.  Shull,  Rev.  Mod.  Phys.,  ^3,  624  ( 1963) . 

4.  H.  Shull  and  F.  Prosser,  J.  Chetn.  Phys.,  40,  233  (  1964) . 

5.  H.  Shull,  j.a.c.s.,  06,  1469  (1964). 

6.  B.  Anex  and  H.  Shull,  "Molecular  Orbitals  in  Chemistry,  Physics,  and 
Biology",  edited  by  P.  0.  Lowdln  and  B.  Pullman,  (Academic  Press,  1964), 

P.  227. 

7.  D.  R.  Hartree,  Proc.  Camb.  Phil.  Soc.,  24,  09  ( I928) ;  V.  Fock,  Zelt.  f. 
Phys Ik,  61,  126  (1930) . 

8.  See,  for  example:  R.  K.  Nesbet,  "Advances  In  Chemical  Physics",  edited 
by  I.  Prlgoglne  ( Interscience  Publisher^  Inc.,  1964),  Vol .  9,  p.  32I, 
and  references  therein. 

9.  P.  0.  Lowdln,  "Advances  In  Chemical  Physics",  edited  by  I.  Prlgoglne, 
(interscience  Publishers,  Inc.,  1959),  Vol.  2,  p.  209. 

10.  H.  Shull,  J.  Chem.  Phys.,  30,  1405  (1959). 

11.  R.  E.  Chrlstof fersen,  J.  Chem.  Phys.,  b\,  960  (1964). 

12.  G.  P.  Barnett,  J.  Llnderberg  and  H.  Shull,  J.  Chem.  Phys.,  4jS,  00  ( 1965) . 

13.  P.  0.  Lowdln  and  H.  Shull,  Phys.  Rev.,  _101,  1750  ( 1956) . 

14.  A.  Macias,  to  be  published. 

15.  A.  Weiss,  Phys.  Rev.,  122,  1026  ( I96I) . 

16.  W.  A.  Lester,  Jr.  and  M.  Krauss,  J.  Chem.  Phys.,  44  ,  207  ( 1966) . 

17.  H.  Preuss,  Zelt.  f.  Natur.,  20,  21  ( 1965) . 

18.  B.  D.  Joshi,  J.  Chem.  Phys.,  44,  3627  ( 1966) . 

19.  H.  Shull  and  P.  0.  Lowdln,  J.  Chem.  Phys.,  2J,  I565  (1955). 

20.  G.  V.  Nataroff  and  J.  0.  Hlrachfelder ,  J.  Chem.  Phys.,  715  ( 1963)  . 


43 


21.  S.  Payer imhoff,  J.  Cham.  Phys.,  42,  996  11965). 

22.  W.  Kolos  and  C.  C.  J.  Roothaan,  Rev.  Hod.  Phys.,  2£»  219  (I960). 

23.  P.  0.  Lowdin,  "Advances  in  Chemical  Physics",  I.  Prigogine,  ed., 
(Interscience  Publishers,  1959),  vol.  2,  p.  207. 

24.  A.  W.  '.Veiss,  Phys.  Rev.,  1826  (1961). 

25.  H.  Conroy,  J.  Chem.  Phys.,  41,  1341  (1964). 

26.  H.  Michels,  J.  Chem.  Phys.,  44,  3834  (1966). 

27.  L.  Wolniewicz,  J.  Cham.  Phys.,  42,  1087  (1965). 

28.  A.  D.  McLean,  A.  <Veiss,  and  M.  Yoshimine,  Rev.  Mod.  Phys.,  32,  211 

(1960).  y  ^ 

29.  This  effect  has  also  been  noticed  (without  comments  concerning  the 
bonding  description)  by  A.  G.  Pearson,  R.  D.  Poshusta  and  J.  C. 
3rowne,  J.  Chem.  Phys.,  44,  1815  (1966). 

30.  G.  S.  Handler  and  J.  R.  Arnold,  J.  Chem.  Phys.,  22»  144  11957). 

31.  M.  Schwartz,  Ph.D.  Thesis,  Vanderbilt  University,  1966. 


44 


II.  MICROWAVE  SPECTROSCOPY  OF  TRANSIENT  ENERGETIC  SPECIES 
1.  INTRODUCTION 

Our  microwave  spectroscopic  studies  have  as  their  goal  the  elucidation 
of  the  structure,  reactivity,  and  energetics  of  gaseous  free  radicals  and 
excited-state  molecules,  wfhile  a  considerable  amount  of  work  has  been 
carried  out  by  optical  absorption  and  emission  studies,  very  little  micro- 
wave  spectroscopy  has  been  performed  on  these  unstable  species.  The 
advantages  of  this  technique  may  be  sunrtarized  as  follows. 

(a)  The  microwave  technique  is  the  most  suitable  for  structural 
determinations.  This  information  is  needed  for  the  calculation  of  the 
thermodynamic  properties  of  the  unstable  species,  and  for  an  understanding 
of  the  mechanism  of  their  formation  and  reactions. 

(b)  The  microwave  technique  is  a  suitable  means  for  detection  of 
low  concentrations  of  gas-phase  species. 

(c)  Microwave  spectroscopy  lends  itself  well  to  double-irradiation 
methods;  for  example,  it  is  possible  to  irradiate  a  sample  with  ultraviolet 
radiation  while  simultaneously  observing  its  microwave  spectrum. 

In  the  following  sections  we  have  summarized  the  microwave  studies 
that  have  been  undertaken.  Much  of  the  work  to  date  has  involved  the  design, 
construction  and  testing  of  a  spectrometer  suitable  for  the  proposed 
studies,  so  these  aspects  have  been  discussed  in  some  detail.  Some 
preliminary  studies  of  biacetyl  photosensitization  have  also  been  dis¬ 
cussed,  and  our  initial  attempts  to  produce  and  detect  free  radicals  have 
been  outlined. 
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2.  SPECTROMETER  CONSTRUCTION 

a.  The  Complete  Spectrometer 

In  Figure  11-1  a  block  diagram  of  the  spectrometer  is  presented.  In 
the  following  we  will  give  more  detail  about  seme  of  the  components.  The 
klystrons*  waveguide  components  and  detector  are  standard  items.  At 
present,  our  frequency  range  is  8  GHz  to  30  GHz.  The  klystron  power 
supply  is  a  PRD  universal  type  812  and  is  capable  of  supplying  any  klystron 
that  we  might  need.  Lock-in  (or  phase  sensitive)  detection  is  accomplished 
by  means  of  a  PAR  model  JB-5  lock-in  amplifier  which  operates  between  10  cps 
and  100  kc/s.  Some  pre-amplification  is  performed  by  means  of  a  narrowband 
amplifier  tuned  at  the  reference  frequency.  Spectrum  observation  is 
performed  by  use  of  either  the  Tektronix  type  502A  dual-beam  oscilloscope 
or  the  Bausch  and  Lomb  VOM-6  recorder.  Signal  modulation  of  paramagnetic 
species  is  produced  by  the  Zeeman  modulator  built  and  designed  in  our 
labs.  This  unit  produces  a  magnetic  field  of  up  to  5  to  10  gauss  in  the 
solenoid  at  a  frequency  of  1  kc/s.  Rise  and  fall  times  of  the  1  kc/s 
square-wave  are  about  and  the  bottom  of  the  square  wave  is  zero-based. 
Figure  11-1  presents  a  circuit  diagram  of  this  unit.  The  sample-cell, 
quartz  rod  transmistion  system  and  modulation  coil  are  discussed  in  more 
detail  in  the  following  section. 

b.  Sample  Cell  Design 

Figure  11-2  presents  a  sketch  of  the  cell  and  quartz  rod 
transmission  system.  The  principles  of  this  type  of  dielectric  rod 
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Figure  11-2.  Sw?!«  cell. 


A  (&C-01 
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transmission  system  have  been  discussed  by  Brackett  et  al.1  By  choosing  the 
rod  diameter  properly,  the  system  can  be  made  to  transmit  useable  microwave 
power  over  any  frequency  range  of  interest.  The  hybrid  TEfc.  mode  is  easily 
excited  from  a  waveguide  system  by  means  of  standard  microwave  horns. 

Figure  11-4  presents  a  graph  which  shows  the  useable  power  (that  which  is 
outside  of  the  quartz  rod  but  within  the  quartz  cell)  as  a  function  of 
frequency  for  three  different  rod  diameters. 

At  the  present  time  the  teflon  windows  are  sealed  to  the  quartz  tube 
by  means  of  Apiezon  W  wax,  while  the  vacuum  seal  between  the  rod  and  the 
window  is  produced  by  paraffin.  The  use  of  0-r,ngs  is  being  contemplated  in 
place  of  Apiezon  IV,  but  at  present  this  modification  has  not  been  undertaken. 

The  Vycor  tube  is  the  major  component  of  the  sample  cell;  its  nominal 
dimensions  are  3  feet  x  l”/«6  inch.  Since  Vycor  transmits  into  the  near 
UV,  external  irradiation  of  samples  is  possible.  The  cell  is  attached  to 
a  high-speed  pumping  system  so  that  continuous  flow  techniques  are  possible. 

Zeeman  modulation  is  produced  by  the  solenoid  wrapped  around  the  sample 

cell.  The  coil  is  made  from  16  gauge  cop per  wire  and  has  about  4  turns/inch 
over  a  length  of  30  inches. 

c.  Spectrometer  Tests 

Before  *ntioning  the  results  using  Zeeman  modulation,  a  fee  comm.nts 
about  Stark  (electric  fi.id  modulation)  seem  in  order.  Using  a  smaller 
diameter  sample  cell,  a  5  Kc/s  square  wave  voltage  was  applied  to  two 
aluminum  plates  <3"  X  30"  X  '/,.■)  placed  parallel  to  the  avis  of  the  cell 
and  to  each  other.  Sith  our  square  wave  generator  it  was  possible  to 
obtain  an  electric  field  of  500  volt/cm.  Some  OCS  absorption  lines  in  the 
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Flfur*  11-4,  OUlaetric  rod  tranwi 
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vicinity  of  24325  Mc/s  were  observed  with  this  arrangement,  but  the 
sensitivity  was  very  poor.  It  was  only  1  X  1CT6  cm'1  (minimum  detectable 
absorption  coefficient),  compared  to  about  5  X  1C)"9  for  standard  Stark 
spectrometers.  This  very  poor  sensitivity  was  evidently  caused  by  poor 
modulation  resulting  from  charges  built  up  on  the  Vycor  surfaces.  It  seems 
evident  that  the  arrangement  is  not  satisfactory  so  Stark  modulation 
methods  will  be  of  little  use.  Because  we  are  interested  in  paramagnetic 
species,  no  further  work  along  this  line  is  contemplated. 

Our  first  satisfactory  Zeeman  modulation  spectrum  has  been  obtained  by 
using  N02  as  the  test  sample.  Six  of  the  eight  hyperfine  components^  of 
the  J  =  232>22  -»  j  =  24Ij23  transition  at  26,600  MHz  have  been  observed 
with  S/N  of  about  50/1.  This  places  our  present  spectrometer  sensitivity 
between  1  X  IQ'9  and  1  X  10-.  These  res  s  are  very  encouraging,  since 
we  now  have  under  construction  a  quality  Zeeman  modulator  which  will 
operate  in  the  10-20  Kc/s  range.  Assuming  that  crystal  noise  predominates, 
the  new  unit  should  produce  a  sansitivity  increase  of  a  factor  of /|0  or 
v'/20>  i.e.,  about  3  or  4.  In  fact,  greater  enhancements  are  expected, 
since  the  rise  and  fall  times  will  be  only  about  3*  f3r  the  new  unit,  whereas 
>e  unit  of  Figure  11-3  has  very  poor  rise  and  fall  characteristics.  Table 
IM  gives  a  summary  of  the  specifications  for  the  Zeeman  modulation  unit 
now  under  construction}  the  first  tests  are  expected  to  be  made  in 
December,  1967.  It  seems  not  unreasonable  to  expect  sensitivities  of 
greater  than  1  X  10"a  cm"1  when  the  new  modulator  is  put  into  operation. 

Experiments  for  the  detection  of  OH  radical  are  now  underway,  since 
these  will  provide  a  more  stringent  test  of  our  apparatus.  The  ^-doublet 
transitions  which  occur  at  about  23822  MHz3  will  be  utilized  for  detection 
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TABLE  1 1 -I .  New  Zeeman  Modulator  Specifications 


Power  required: 

Maximum  coil  inductance: 
Frequency: 

Rise  and  fall  time: 

Dr  ive  current: 

'.'/aveform  flatness: 
Zero-basing  of  waveform: 


75  watts  max 
4‘>0  pH 

Variable  to  20  KHz 
i  a  sec  each  at  10  KHz 
4  amps  max 
Better  than  1  y% 

To  within  1  ma 
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purposes.  OH  radicals  will  be  produced  by  use  of  a  2500  MHz  microwave 
discharge  by  one  of  the  following  systems: 


discharge 

h2  ♦  o2 

- ^  >  OH  +  ... 

(D3 

H  +  N02 

- >  OH  «•  NO 

(2)4 

HaC  ♦  H2G2 

discharr^  UH  +  _ 

(3)5 

The  radical  production  will  be  made  continuous  by  using  the  rapid  flow 
technique  previously  mentioned. 

'"e  plen  to  detect  SO  radicals  with  our  Zeeman  spectrometer  also; 
previous  workers4  have  observed  several  transitions  using  Stark  modulation, 
but  no  magnetic  studies  were  performed. 

3.  INITIAL  STUDIES 

a.  free  Radicals 

The  first  attempt  at  detecting  a  new  free  radical  will  be  performed  on 
the  Cf3  species.  CF^CA)  has  been  detected  by  Powell  and  Lide  as  a 
product  of  electrodeless  discharges  in  a  variety  of  fluorinatcd  molecules 
including  (CF3)2C0.  '.7e  will  search  for  the  symmetric  rotor  spectrum  of  Cf3 

produced  by  flowing  perf luoroacetone  through  a  microwave  discharge.  It 
seems  highly  unlikely  that  this  free  radical  will  have  zero  dipole  moment, 
and  it  is  expected  to  exhibit  its  J  =  0—>  J  =  1  transition  in  the 
neighborhood  of  20  GHz.  It  should  be  noted  that  absorption  lines  of  non- 
paramagnetic  species  will  not  be  observed  with  the  Zeeman  system,  so  that 
the  parent  molecule  spectrum  will  cause  no  difficulty. 
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A  search  will  also  be  carried  out  for  the  spectrum  of  the  K)  (2n) 

n 

radical.  This  species  was  recently  observed  in  emission  by  passing 
phosphorus  vapor  and  oxygen  through  a  2450  MHz  microwave  discharge.  Although 
the  J  =  0 — >  J  =  1  transition  is  beyond  our  frequency  range  (at  about  45  GHz), 
there  is  some  hope  of  observing  magnetic  dipole  transitions  at  lower 
frequencies. 

For  radical  studies  in  which  it  is  desired  to  mix  the  products  of  the 
discharge  with  another  species,  we  use  the  arrangement  shown  in  Figure  11-5. 

This  arrangement  is  useful,  for  example,  in  producing  OH  by  mixing  N0k  with 
H  atoms  produced  by  dissociation  of  H2  in  the  discharge. 

b.  Electronically  excited  Species 

Preliminary  experiments  that  we  have  performed  show  tbit  the  production 

of  triplet  biacetyl  may  e  efficient  enough  to  permit  the  observation  of 

its  microwave  spectrum.  T  ,e  method  is  as  follows.  Triplet  naphthalene  is 

o 

formed  by  relatively  strong  rradiation  at  2537  A  of  ground  state  naphthalene, 
followed  by  intersystem  errs  'ng  from  the  excited  singlet  system.  Collisions 
of  the  second  kind  between  sintiot  biacetyl  and  triplet  naphthalene  produce 
triplet  biacetyl. 

The  above  photosensitization  action  has  been  observed  semi-quant itat ively 

in  the  following  manner.  A  mixture  of  biacetyl  and  naphthalene  has  been 

o 

irradiated  at  2537  A.  The  emission  of  the  sample  has  been  observed  with  a 
Bausch  and  Lomb  monochromator  and  has  been  found  to  contain  only  the  green 
biacetyl  phosphorescence;  no  biacetyl  fluorescence  is  observed.  This 
clearly  indicates  that  energy  transfer  occurs  as  described  above,  and  not 
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Figwra  11-5.  Dlacharga  and  aUing  t»*a. 
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by  transfer  of  excited  singlet  energy  between  the  molecules.  Experiments  to 
determine  the  quantum  efficiency  of  this  process  are  underway  and  will  be 
ready  for  publication  in  a  short  time. 

The  microwave  experiments  will  be  performed  by  passing  a  mixture  of 

biacetyl  and  naphthalene  through  the  cell  while  irradiating  with  mercury 

resonance  lamps.  Crude  calculations  indicate  that  about  1020  photons/sec 
0 

of  2537  A  radiation  will  be  sufficient  to  produce  detectable  amounts  of 
triplet  biacetyl.  Lamps  are  available  which  supply  this  much  radiation,  but 
they  will  have  to  be  specially  fabricated  to  achieve  such  intensities  over 
our  geometric  configuration. 
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III.  RADI  AT  I  ON- INDUCED  CHARGE  STORAGE  IN  ORGANIC  POLYMERS 
1.  INTRODUCTION 

The  primary  interaction  of  gamma  rays  of  energies  between  0.3  and  3  Mev 
with  matter  is  the  Compton  process.  Electrons  scattered  by  this  mechanism 
leave  positive  holes  and  move  through  the  medium  causing  additional  ionization 
and  excitation  until  they  either  recombine  or  are  trapped.  Likewise 
positive  holes  are  free  to  move  in  some  materials  until  they  are  neutralized 
or  trapped.  Charge  carriers  which  are  stabilized  in  a  solid  medium  can  be 
described  as  trapped  in  a  potential  well.  The  charges  may  be  det rapped  by 
thermal  or  photo  excitation.  The  detrapped,  mobile  charges  can  be  used  to 
probe  the  cross-sect  ion  for  retrapping  which  relates  to  the  nature  of  the 
trapping  site,  the  spatial  distribution  of  trapping  sites,  the  mobility 
mechanism  and  the  limiting  density  of  trapped  charges. 

It  is  of  importance  to  study  both  trapped  and  mobile  charges  in  order 
to  properly  assess  the  feasibility  of  storing  large  amounts  of  charge  for 
possible  rocket  propellants.  In  the  present  program  this  is  being  achieved  by 
thermally  generated  currents  in  irradiated  organic  solids.  The  development 
of  electron  paramagnetic  resonance  studies  of  trapped  charges,  reactivity 
studies  of  mobile  ions  in  the  liquid  phase  and  mobility  measurements  of 
electrons  and  holes  in  irradiated  solids  is  also  being  implemented. 

Thermal  heating  can  be  used  to  detrap  charges  in  irradiated  solids.  If 
a  temperature  gradient  is  applied  across  the  sample  the  charges  will 
preferentially  drift  in  one  direction  and  create  an  observable  current. 

Such  currents  are  called  thermoelectric  currents.  Hardtke^  reports  the 
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observation  of  thermoelectric  currents  of  a  few  picoamperes  when  samples  of 
cabal  glass  are  irradiated  by  cobalt-60  and  then  heated  in  temperature  gradi¬ 
ents  approximating  300C/mm.  He  also  reports  observing  the  same  phenomenon 
in  certain  other  dielectrics  including  polystyrene  and  nylon.  Gross2 

measured  similar  currents  under  similar  temperature  gradients  in  borosilicate 
glass. 

Experimentally,  one  can  consider  two  classes  of  organic  systems  to  be 
studied.  In  the  first  class  are  rather  simple  o:yanics  such  as  3-methylpentane, 
triethylamine,  etc.,  which  are  liquids  or  gases  at  room  temperature,  and 
which  must  be  frozen  and  irradiated  at  temperatures  down  to  4°«  in  order  to 
trap  ions.  In  this  class  exploratory  work  is  beginning  on  triethylamine.  The 
second  class  are  organic  solids  which  can  trap  charge  at  room  temperature 
and  above;  these  are  usually  polymeric  systems. 

Charge  storage  in  polymethylmethacrylate  (PmMA)  has  been  investigated 
3 

by  Barnes.  He  found  that  a  voltage  was  generated  across  an  irradiated 
Lucite  sample  when  heated.  This  suggested  that  charges  were  trapped  in  the 
irradiated  Lucite  and  that  this  material  would  be  a  fruitful  system  for 
further  study.  The  results  obtained  here  by  Barnes  in  a  subsequent, 
independent  investigation  bear  this  out.  By  using  pure  PMMA  as  well  as 
PMMA  doped  with  hole  and  electron  traps  we  have  discovered  that  added  electron 
traps  but  not  hole  traps  lead  to  charge  stabilization.  tVe  have  further 
found  that  polymers  containing  built-in  electron  traps,  such  as  polyvinyl 
chloride  (PVC),  stabilize  even  larger  amounts  of  charge.  It  is  of  con¬ 
siderable  significance  that  the  stored  charge  can  be  obtained  as  a 
thermoelectric  current  by  application  of  extremely  small  temperature  gradients 
of  -5°/ cm.  The  requirements  for  charge  stabilization  in  irradiated  organic 
polymers  at  relatively  high  temperatures  are  beginning  to  emerge. 
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2.  EXPERIMENTAL 


a.  Sample  Preparation 


oarnples  of  the  organic  polymers  under  investigation  were  machined  into 
disks,  normally  2.5  cm  in  diameter  and  1  cm  thick,  from  cylindrical  stock. 
This  geometrical  shape  facilitated  the  attachment  of  electrodes  and  the 
maintenance  of  a  clean  insulating  surface  between  the  electrodes. 

The  flat  surfaces  of  the  disks  were  coated  with  a  low  resistance 
conductive  cement  for  plastics,  Eccobond  Solder  56C,  to  serve  as  electrodes. 
Eccobond  Solder  (Emrnerson  and  Cummings,  Inc.)  is  a  silver-colored  paste  to 
which  a  No.  11  catalyst  is  added  for  hardening.  Cure  of  the  coated  surfaces 
was  accomplished  in  4  hours  at  65°C. 

The  Eccobond  solder  also  provided  a  method  of  attaching  3/4  inch  metal 
tabs  to  the  flat  surfaces  of  the  disks.  These  tabs  were  opposite  each 
other  and  allowed  to  protrude  1/2  inch  from  the  edge  of  each  electrode. 

Each  tab  provided  a  means  of  handling  the  sample  by  one  electrode  at  a  time 
without  shortcircuiting  the  two  electrodes. 

The  primary  organic  polymer  studied  was  polymethylmethacrylate  (PIWUIA). 


-  CH2  -  C  - 

i 

C0QCH3 


monomer  to.  IV.  =  100,11 


This  polymer  is  available  commercially  under  the  trade  names  of  Lucite, 
Plexiglas,  Diakon, Perspex  and  Crystalite.  Usually  the  commercial  polymeri¬ 
zation  involves  the  use  of  benzoyl  peroxide  as  a  chemical  initiator  and  some 
benzoyl  peroxide  is  present  in  trace  quantities  in  the  final  product. 
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In  addition  to  commercial  P*/LA,  studies  were  made  of  material  polymerized 
from  pure  monomer.  Various  chemical  compounds  including  carbon  tetrachloride, 
benzene,  chlorobenzene,  benzoyl  peroxide  and  2-methylpentene-1  were  added 
to  the  monomer  in  small  concentrations  before  polymerization.  Carbon 
tetrachloride  is  a  good  electron  scavenger  and  2-methylpentene-1  is  a 
positive  hole  trap.^ 

Other  materials  stud.ed  include  polyvinyl  chloride  (polychloroethylene) , 
polyethylene,  polytatrsr luoroethylene  (Teflon),  and  triethyleneglycol 
dimethacrylate  (Tedma).  All  of  these  materials  were  obtained  from  corrmercial 
stock  except  Tedma  which  was  thermally  polymerized  from  the  pure  monomer. 

It  was  necessary  to  heat  the  polyvinyl  chloride  samples  to  90°C  for  30  minutes 
in  short  circuit  to  remove  residual  charge  carriers  trapped  in  the  material 
by  the  polymerization  process.  This  preheat  treatment  was  not  necessary 
for  any  of  the  other  materials  studied. 

b.  Current  Measurement 

The  disk-shaped  samples  were  held  during  heating  and  current  measurement 
in  a  three-sided  Teflon  sample  holder  (Figure  1 1 1-1)  between  two  brass 
screws,  one  of  which  was  moveable.  Small  copper  disks  were  soldered  to  the 
heads  of  each  brass  screw  to  provide  a  smooth  flat  surface  for  physical  and 
electrical  contact  with  the  sample  electrodes.  Banana  plugs  were  used  to 
connect  the  brass  screws  with  the  remainder  of  the  current  measuring  circuit. 

An  additional  electrical  connection  was  made  to  the  metal  tab  on  each  sample 

electrode  by  a  copper  wire  soldered  to  the  brass  screw  and  spots  soldered 
to  the  tab. 
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FIgur.  1 1 1-1.  T«flon  holder. 


During  measurements,  the  sample  holder  rested  on  a  brass  platform  with 
on.  electrode  connected  to  the  preamp  input  plug.  The  other  electrode  was 
returned  by  a  heavy  multistranded  wire  to  the  electrometer  ground.  After 
placement  of  the  sample  holder  in  the  oven  a  grounded  brass  cylinder  was 
moved  over  it  to  serve  as  a  guard  ring  and  to  control  the  flow  of  heated  air 
circulated  by  the  oven  fan.  The  brass  cylinder  can  be  considered  to  be  a 
Faraday  cage  and  therefore  supplies  capacitive  and  inductive  shielding  to 
the  sample  and  sample  holder. 

The  oven  fan  v*s  centered  in  the  top  of  the  oven  directly  over  one  of  the 
two  toroidal  shaped  800  watt  oven  heaters.  The  other  oven  heater  was 
located  in  the  bottom  of  the  oven.  The  hot  air  circulated  by  the  fan  moved 
downward  heating  the  grounded  brass  screw  which  protruded  slightly  from  the 
end  of  the  brass  cylinder.  The  other  brass  screw  was  connected  to  the 
input  plug  which  acted  as  a  heat  sink  thus  cooling  that  side  of  the  sample 
holder  and  sample.  This  heating  and  cooling  resulted  in  the  temperature 
gradient  utilized  in  this  investigation. 

The  input  plug  consisted  of  a  12  inch  long,  1  inch  diameter  brass  pipe, 
three  teflon  insulators  and  a  heavy  copper  wire  center  conductor.  The 
purpose  of  this  plug  was  to  prevent  heating  of  the  electrometer  preamplifier. 
The  outer  brass  cylinder  was  connected  to  the  electrometer  ground  and  helped 
reduce  noise  in  the  system. 

c.  Temperature  Measurement 

Figure  MI-2  is  a  block  diagram  of  the  thermoelectric  current  measuring 
circuit.  The  rate  of  heating  in  the  oven  was  controlled  by  a  Varian 
Aerograph  Model  326  linear  temperature  progranwer.  Linear  heating  rates  from 
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Figure  l 11-2.  Circuit  block  diacr, 
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2°C/min  to  40°C/min  can  be  selected.  For  most  measurements  a  rate  of  6°C/min 
was  employed. 

The  temperature  of  the  oven  was  sensed  by  a  copper-constantan  thermo¬ 
couple  soldered  to  the  small  copper  disk  which  made  contact  with  the  hotter 
sample  electrode.  This  was  on  the  grounded  stationary  electrode  of  the 
sample  holder.  The  leads  of  this  thermocouple  were  spliced  with  a  plug- jack 
connection  inside  the  oven  so  the  sample  holder  was  removable.  The  potential 
difference  generated  at  the  thermocouple  junction  was  measured  by  a 
Digitec  Model  564  digital  read-out  thermometer.  This  provided  an  accurate 
continuous  measurement  of  the  temperature  of  the  grounded  sample  electrode. 

The  signal  from  this  thermocouple  was  also  used  as  the  input  for  the  X-axis 
of  the  X-Y  recorder. 

The  thermoelectric  currents  which  were  obtained  by  heating  the  sanples 
of  organic  polymers  were  measured  with  a  Cary  Model  31  vibrating-reed 
electrometer .  An  input  resistance  of  1010  ohms  was  used  for  all  measure¬ 
ments  and  a  signal  proportional  to  the  magnitude  of  the  current  was 
generated.  This  signal  was  plotted  against  temperature  on  a  MoseleyModel 
2  FRAM  X-Y  recorder.  Simultaneously*  plots  of  current  versus  time  and 
temperature  versus  time  were  made  on  a  Moseley  Model  7100A  two  pen,  time 
base  recorder.  Integration  of  the  current  versus  time  plot  yielded  values 
proportional  to  the  normalized  total  charge  released  in  the  sample. 

At  hifsh  temperatures  thermal  stresses  in  the  samples  produce  noise  which 
mask  the  thermoelectric  current.  This  noise  occurs  at  a  high  frequency 
and  can  be  filtered  out  to  leave  only  the  slowly  changing  thermoelectric 
current.  Figure  IM-3  shows  the  filtering  circuit  used  on  the  recorder  input 
signal.  It  is  a  simple  RC  circuit  with  a  variable  time  constant,  0-10  seconds. 
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Figure  1 1 1 —3 •  Noise  filtering  circuit. 
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For  a  very  fast  rising  peak  a  slight  distortion  will  result  but,  assuming  the 
peak  is  symmetrical,  the  area  under  the  curve  will  be  unchanged. 

d.  Temperature  Gradients 

The  temperature  gradient  across  the  sample  was  measured  by  using 
Eccobond  solder  to  attach  copper-constantan  thermocouples  directly  to  the 
conducting  electrodes  of  an  unirradiated  sample.  The  potential  generated 
by  each  thermocouple  was  selected  by  a  switching'  arrangement  and  measured 
by  a  Leeds  and  I'lorthrup  potentiometer.  Standard  conversion  tables  for 
thermocouples  were  used  to  correlate  the  observed  potentials  with  the  corres¬ 
ponding  temperatures.  Values  for  the  two  electrodes  were  subtracted, 
divided  by  the  sample  thickness,  and  than  plotted  against  time  to  obtain 
temperature  gradient  curves.  Figure  111-4  shows  the  results  of  such  a 
temperature  gradient  measurement  on  a  sample  1  cm  thick.  Figure  1 1 1-5  shows 
a  similar  temperature  gradient  measurement  on  a  thin  sample  0.1  cm  thick. 

In  order  to  test  the  postulate  that  the  temperature  gradient  caused 
the  flow  of  measured  thermoelectric  current,  measurements  were  made  on 
samples  heated  in  the  absence  of  a  temperature  gradient.  Some  experimentation 
showed  that  by  moving  the  sample,  sample  holder  and  brass  guard  cylinder  to 
the  center  of  the  oven  a  position  could  be  determined  at  which  both  sides  of 
the  sample  holder  were  heated  at  equal  rates.  Thus  copper  wire  extensions 
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Figure  1 1 1-5.  Temperature  gradient  measurement  on  Lucite  0.1  cm  thick. 
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to  the  input  plug  and  electrometer  ground  were  necessary  to  prevent  any 
appreciable  conduction  of  heat  away  from  the  sample  holder.  A  plot  of  the 
t'jmperature  difference  of  the  sample  electrodes  when  heated  in  this 
configuration  is  shown  in  Figure  I  I  1-6. 

a.  Polymer i zat ion  Techniques 

Inhibitor  free  methylmethacrylate  (MMA)  monomer,  99.95#  pure,  was 
obtained  from  the  Rohm  and  Haas  Company.  Polymerization  of  this  monomer 
occurs  by  a  free  radical  process.  Unreacted  monomer  molecules  ado  to 
monomer  free  radicals  to  form  the  polymer  chain.  The  process  is  known  as 
addition  polymerization.^  If  the  monomer  free  radicals  are  the  result  of 
the  presence  of  a  chemical  initiator  (I)  the  process  can  be  represented  by 
the  following  reactions: 


I 


2R' 


R’  +  M  - M/ 

V  ♦  M  - >  M2* 


M  * 

x 


Vi 


The  chain  continues  to  propagate  until  the  monomer  molecules  become  depleted, 
the  terminal  free  radicals  combine  with  other  species,  or  the  viscosity 
increases  until  monomer  molecules  can  no  longer  move  to  the  end  of  the  chains. 

Polymerization  of  monomeric  MMA  can  be  initiated  by  ultraviolet  light, 
traces  of  c '.ygen  or  peroxides,  elevation  of  temperature,  or  by  the  addition 
ct  certain  chemical  initiators.  Polymerization  with  ultraviolet  light 

(photopolymerization)  was  accomplished  in  a  Rayonet  U.  V.  Reactor.  The  light 

.  o 

source  consists  of  sixteen  3500  A  mercury  vapor  fluorescent  lanps  arranged 
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Figure  I  I  1-6.  Zero  temperature  gradient  measurement. 
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in  a  cylindrical  configuration  with  a  combined  intensity  given  by  the 
manufacturer  of  1.5  to  5  X  1016  photons/sec/cm3.  The  ultraviolet  irradiation 
was  carried  out  in  cylindrical  pyrex  glass  tubes.  A  24-hour  exposure  was 
necessary  to  obtain  a  hard  polymer. 

Thermal  polymerization  of  MMA  required  4 2  hours  at  a  temperature  of  90°C. 

7 

Riddle  recommends  a  variation  of  temperature  in  a  cycle  ranging  from  40  to 
110  C  in  order  to  control  the  polymerization  process  during  its  accelerated 
stages.  This  method  is  now  being  explored. 

g 

The  effect  of  oxygen  on  the  polymerization  of  MMA  has  been  investigated 
and  found  to  have  a  definite  inhibiting  effect  on  the  polymerization.  It  has 

been  shown  that  a  low  molecular  weight  polymeric  peroxide  is  formed  with  MMA 

.  9 

which  has  the  following  probable  structure. 

CH, 

i 

-ch2-c-o-o- 

i 

ccgch3 

J  n 

polyperoxide 

This  polyperoxide  then  decomposes  to  give  formaldehyde  and  methylpyruvate.  In 

addition  the  oxygen  probably  acts  as  a  chain  terminator  in  the  early  stages 

of  polymerization.  The  polymeric  peroxide  can  initiate  polymerization  which 

7 

may  sometimes  be  extremely  rapid. 

Oxygen  was  removed  from  the  monomer  prior  to  polymerization  by  the 
freeze  -  pump  -  thaw  method  for  all  samples  reported  herein.  It  has  been 
found  that  the  monomer  can  be  equally  effectively  deoxygenated  by  purging 
for  15  minutes  with  dry  nitrogen.  All  future  samples  will  be  deoxygenated  by 
purging  since  it  is  more  rapid  than  the  freeze  -  pump  -  thaw  method. 
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As  mentioned  previously  benzoyl  peroxide  is  the  most  commonly  used 

chemical  initiator  in  the  commercial  polymerization  of  MMA.  A  tertiary  amine 

act i vator ,  such  as  N,M-dimethylani  1  ine,  N,N-dimethyl-p-toluidine  or 

tr i hexyl amine  is  sometimes  used  with  benzoyl  peroxide  to  speed  up  the  reaction 
7 

at  room  temperature. 

All  data  on  synthesized  polymers  presented  in  this  report  were  obtained 
from  photopolymer i zed  samples  unless  the  material  added  was  benzoyl  peroxide 
which  itself  served  to  catalyze  the  polymerization.  Thermal  polymerization 
methods  arepresently  being  explored. 

f.  Sample  Irradiation 

All  samples  were  irradiated  in  a  2900  curie  U.  3.  Nuclear  GR-9  Co-60 
Irradiator  at  a  dose  rate  of  approximately  .5  Mrad/hr.  The  source  is  arranged 
in  a  cylindrical  configuration  around  the  4"  diameter  sample  chamber  to  pro¬ 
vide  approximately  isotropic  irradiation.  Samples  were  supported  by  one  of 
the  metal  tabs  attached  to  the  sample  electrodes,  and  were  positioned  in  the 
center  of  the  irradiation  chamber  to  assure  uniform  exposure.  The  ambient 
temperature  of  the  source  is  about  35°C. 
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3.  RESULTS 

a.  Currant  Measurement 

A  typical  current  vs.  time  plot  of  a  1  cm  thick  Lucite  sample  is  shown 
in  Figure  II 1-7.  The  sample  was  heated  linearly  at  a  rate  of  6°C/min 
after  a  dose  of  0,5  Mr ad.  The  curve  is  approximately  symmetrical  with  a 
peak  of  5.9  picoamps  (pa)  at  65°C.  Figure  II 1-6  is  a  similar  plot  for  a 
PVC  sample  measured  under  the  same  conditions.  The  peak  current  is  31  pa 
and  occurs  at  a  temperature  of  86°C.  A  second  heating  of  each  sample 
resulted  in  no  measurable  current. 

All  temperatures  mentioned  were  measured  at  the  higher  temperature 
sample  electrode.  For  temperatures  above  90°C  in  Lucite  electronic  noise 
was  observed  on  the  recorder.  This  temperature  corresponds  to  the  glass 
transition  temperature  (Tg)  for  PtoMA.  The  filtering  circuit  mentioned 
previously  allowed  the  measurement  of  the  thermoelectric  current  in  the 
presence  of  noise  to  a  temperature  of  120°C.  Heating  samples  above  150°C 
caused  gas  bubbles  to  form  in  the  material  and  caused  the  samples  to  swell. 

Integration  of  the  current  vs.  time  curves  does  not  yield  values  equal 
to  the  actual  total  charge  released  in  the  sample  but  rather  gives  values 
proportional  to  the  total  charge  release.  This  is  because  each  charge 
carrier  released  moves  only  a  fraction  of  the  distance  between  the 
electrodes  before  it  is  neutralized,  and  therefore  contributes  only  a 
fraction  of  a  charge  carrier  to  the  total  current  measured  in  the  external 
circuit.  Integration  of  all  curves  to  obtain  normalized  values  propor¬ 
tional  to  the  total  charge  released  in  a  sample  were  carried  out  to  a 
temperature  of  120°C. 
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Figure  111-8.  Thermoelectric  current  plot  for  PVC. 

The  asymmetry  is  often  observed  and 
its  origin  is  unknown. 
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b.  Dose  Effects 
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Samples  of  Lucite  1  cm  thick  were  given  doses  from  0.12  Mrad  to  20 
Mrad.  These  samples  were  then  heated  linearly  at  a  rate  of  6°C/min  while 
the  currents  were  measured.  Results  of  the  measurements  are  shown  in 
Table  1 1 1-1  and  Figure  1 1 1-9.  The  total  charge  released  from  the  sample 
increases  with  dose  to  about  1  Mrad.  It  then  appears  to  remain  relatively 
constant  to  at  least  a  dose  of  20  Mrad.  This  constant  value  of  normalized 
total  charge  released  suggests  a  saturation  of  trapped  charge  carriers  in 
the  material  at  quite  low  doses. 

c.  Effect  of  Heating  Rate 

A  change  in  the  rate  of  heating  a  sample  primarily  changes  the  tempera¬ 
ture  gradient  imposed  on  the  sample.  The  temperature  gradient  influences 
the  movement  of  charge  carriers  predominantly  in  one  direction  and  therefore 
the  recombination  of  charge  carriers  of  opposite  sign.  Table  lll-ll 
contains  the  important  parameters  of  thermoelectric  current  peaks  measured 
at  heating  rates  from  2°C/-nin  to  12°C/min.  The  variation  in  normalized 
total  charge  with  heating  rate  is  shown  graphically  in  Figure  111-10  while 
Figure  111-11  illustrates  the  relation  of  heating  rate  to  peak  current. 

Both  parameters  increase  linearly  with  heating  rate  although  the  effect  on 

peak  current  was  much  greater.  All  samples  were  1  cm  thick  and  were  given 
a  dose  of  0.5  Mrad. 

d.  Zero  Temperature  Gradient 

A  position  in  the  oven  had  been  determined  as  described  previously 
(Figure  1 1 1-6)  at  which  a  sample  could  be  heated  with  a  very  small  temperature 
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Table  III -I:  Effect  of  Absorbed  Dose  on  Thermal 
Induced  Charge  Release  In  Luclte 


Heating  Rate:6°C/min. 

Samples  1 

cm  Thick 

Exposure 

Total 

Peak 

Temperature 

Normalized 

Time 

Dose 

Current 

at  Peak 

Total  Charge* 

Hours 

Mrad 

picoamps. 

°C 

arbitrary  units 

0.25 

0.12 

3-9 

65 

119.9 

0.5 

0.25 

4.8 

65 

127.5 

1 

0.5 

5-9 

65 

168.0 

2 

1.0 

6.2 

62 

155.6 

10.25b 

5.1 

5.8 

66 

158.1 

O.45 
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20b 

10 

5.7 

66 

I6I.3 

O.65 

120 

4ob 

20 

5.6 

62 

158.1 

1.3 
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Figuro  I II -9.  Thormal  induced  charg®  releaa®  vs.  absorbed  dosa. 


Table  III -II:  Effect  of  Heating  Rate  on  Thermal 
Induced  Charge  Release  in  Lucite 


Absorbed  Dose: 0.5  Mrad 

Samples  1 

cm  Thick 

Heating 

Peak 

Temperature 

Normalized 

Rate 

Current 

at  Peak 

Total  Charge* 

sC/min 

plcoamps 

°C 

arbitrary  units 

2 

1.9 

53 

108.0 

4 

2.8 

59 

100.5 

6 

5.9 

65 

168.0 

8 

6.5 

65 

141.6 

10 

7.6 

69 

155-2 

12 

10.4 

68 

193.8 

* 


a 


Curves  integrated  to  120 °C 
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Figure  111-10.  Effect  of  heatir^  rate 


on  noraalf serf  total  charge. 


2  4  6  8  10  12 

HEATING  RATE  CC/MIN.) 

Figure  111-11.  Effect  of  heating  rate  on  peak  current. 
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difference  between  the  sample  electrodes.  In  order  to  determine  the 
importance  of  the  temperature  gradient  on  the  current  measurements,  an 
experiment  was  performed  with  a  sample  in  this  position.  The  result  of 
this  heating  is  shown  in  Figure  111-12.  The  Lucite  sample  was  1  cm  thick 
and  had  received  0.5  Mrad  of  gamma  irradiation.  The  oven  was  heated 
linearly  at  a  rate  of  6°C/min.  Both  the  peak  height  and  area  under  the 
curve  were  reduced  by  factors  of  approximately  10  and  20, respectively, 
compared  with  samples  measured  under  similar  conditions  but  in  the  normal 
sample  position.  This  is  evidence  that  a  temperature  gradient  is  essential 
to  the  measurement  of  the  observed  thermoelectric  currents. 

e.  Successive  Irradiations  on  the  Same  Sample 

All  measurements  made  in  this  investigation  were  performed  on  samples 
with  no  known  previous  irradiation  history  with  the  exception  of  those 
reported  in  this  section.  Table  Ill-Ill  gives  the  parameters  for  two 
different  Lucite  samples  irradiated  and  measured  three  successive  times. 
Similar  measurements  were  made  on  two  PVC  samples.  Table  1 1 1  —  I V .  The 
radiation  was  expected  to  produce  products  which  might  act  as  additional  or 
deeper  trapping  sites.  An  increase  in  normalized  total  charge  is  apparent  in 
most  samples. 

f.  Variations  of  Sample  Thickness 

The  results  of  measuring  the  thermoelectric  current  in  samples  of 
different  thickness  ara  shown  in  Figures  111-13  and  111-14*  In  Figure  111-13 
the  current  vs  temperature  plots  for  selected  samples  between  1.0  and 
0.17  cm  thick  are  superimposed  and  the  same  curves  for  samples  between  0.10 
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Figure  111-12.  Thermoelectric  current  plot  for 
zero  temperature  gradient. 
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Table  111-111:  Successive  Irradiations  and 
Measurements  of  Thermoelectric  Currents  In  Lucite 


Absorbed  Doae:0.5  Mrad  Heating  Rate: 6°C/mln. 


Trial 

no. 

Peak 

Current 

Temperature 
of  Peak 

Normalized 
Total  Charge 

picoamps 

8C 

arbitrary  units 

Sample  no.  229 

1 

5.5 

85 

178.9 

2 

5.7 

86 

180.9 

3 

Sample  no.  230 

5.5 

84 

192.1 

1 

4.8 

86 

160.8 

2 

5.7 

86 

180.9 

3 

5.6 

83 

210.5 

86 


Table  III -IV:  Successive  Irradiation  and  Measurement  of 
Thermoelectric  Currents  in  Polyvinyl  Chlorit’e 


Absorbed  Do,e;0.5  Mrsd  Heating  Rate:6»C/min. 


Trial 

no. 

Peak 

Current 

Temperature 
of  Peak 

Normalised 
Total  Charge 

picoampo 

°C 

arbitrary  unita 

Sample  no.  23I 

I 

1 

31 

85 

67.1* 

2 

32 

85 

7k.6 

3 

38 

86 

76.3 

1* 

void 

void 

void 

5 

1*7 

82 

87.8 

6 

1*8 

82 

85.U 

Sample  232 

1 

31 

86 

67.5 

2 

32 

81* 

7^.5 

3 

35 

86 

81.5 

CURRENT 


TEMPERATURE  ( • c  ) 


CURRENT  fPICO  AMPS) 


ri*'jr®  ll,_U*  Current  profile  for  aempl*  tMckn*s.«« 
from  0.10  cm  to  0.02  cm. 
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and  0.02  cm  thick  ara  illustrated  in  Figure  lll-U.  All  measurements  were 
made  after  a  dose  of  0.5  Mr ad  and  for  a  linear  heating  rate  of  6°C/min. 
Integration  of  the  current  vs  temperature  plots  was  performed  to  obtain 
relative  values  proportional  to  the  net  charge  released  in  the  sample.  The 
net  area  ms  obtained  by  subtracting  the  area  under  the  positive  and 
negative  current  peaks.  These  values  were  plotted  against  sample  thickness 
in  Figure  111-15.  Examination  of  this  figure  indicates  a  maximum  net  charge 
release  for  a  sample  thickness  of  approximately  0.4  cm  and  a  zero  net 
charge  release  for  a  thickness  of  about  0.09  cm. 

g.  Decay  Studies 

A  study  was  made  of  the  decay  of  the  total  charge  released  in  Lucite 
at  various  temperatures.  Samples  1.0  cm  thick  were  given  doses  of  0.5  Mrad 
and  after  decay  at  various  constant  temperatures,  they  were  measured  at  a 
linear  heating  rate  of  6°C/min.  Table  lll-V  shows  the  results  of  these 
measurements. 

The  normalized  total  charge  vs  decay  time  for  Lucite  samples  is  shown 
in  semilog  plots  in  Figure  111-16  for  selected  temperatures  from  0°C  to 
50°C.  All  the  decay  curves  appear  to  contain  two  components,  a  rapid  initial 
decay  and  a  slower  long  term  decay. 

Decay  constants  were  obtained  from  the  slope  of  linaar  portions  of  the 
Lucite  decay  curves  and  were  plotted  on  semi  log  paper  against  the  inverse 
of  the  absolute  temperature  to  obtain  an  Arrhenius  plot,  Figure  111-17. 

The  slope  of  this  line  is  equal  to  the  activation  energy  (Ea)  divided  by  the 
universal  gas  constant  (R).  For  Lucite  the  activation  energy  was  calculated 
to  be  20.7  kcal/mole  (0,9  ev).  Preliminary  decay  measurements  on  PVC  show 
that  trapped  charges  are  much  more  stable  in  it.  Detailed  decay  studies  on 
PVC  are  in  progress. 
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Tabic  111 -Vs  Decay  of  Trapped  Charge  Carriers  in  Lucite 
at  Various  Tanperstures' 


Samples  1  cm  Thick 

Absorbed  Dose: 0.5  Mrad 


Decay 

Tine 


Peak 

Current 


perature  o“C 
5 


Heating  Rate  6*0 /min 


Temperature 
of  Peak 


Normalized 
Total  Charge 


arbitrary  units 


TOTAL  CHANGE  (orb.  unit*) 


DECAY  TIME  (hoar*) 

Figure  111-16.  Decay  of  thermoelectric  currente  in  Luclte  eenplee. 
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Figure  111-17.  Arrhenius  plot  for  thermoelectr ic  current  decay. 
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h.  Effect  of  Additives 

Pure  PMMA  was  photopolymerized  with  and  without  various  selected 
additives  as  described  previously.  Samples  were  prepared  from  these 
synthesized  materials.  The  samples  were  given  a  dose  of  0.5  Mrad  and 
subsequently  were  heated  at  a  linear  rate  of  6°C/min  for  current  measure- 
ment.  The  results  of  this  study  are  given  in  Taole  lll-VI.  Pure  PIUMA 
with  no  additives  produces  a  very  snail  thermoelectric  current  while  the 
presence  of  small  concentrations  of  electron  traps  like  carbon  tetrachloride 
or  benzoyl  peroxide  increase  this  current  by  a  factor  of  at  least  10. 

Similar  concentrations  of  benzene  and  chlorobenzene  do  not  increase  the 
current  over  that  in  pure  PMUA.  The  hole  trap>  2-methylpentene— 1 »  causes 
a  twofold  increase  in  the  current  for  one  run;  further  experiments  are 
necessary  to  ascertain  if  this  increase  is  significant.  The  relatively 
large  current  from  PVC  samples  is  not  unexpected  in  view  of  the  results  of 
CCI4  as  an  additive. 

A  measurement  on  irradiated  polyethylene  resulted  in  a  negligible 
thermoelectric  current.  Samples  of  Teflon  and  triethyleneglycoldimeth- 
acrylate  (Tedma)  also  produced  nearly  zero  thermoelectric  currents  after 
irradiation  and  heating. 
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Table  III -VI: 


Material 


neu 

Jffcpl,  i 

A^La  J 


AMA 


Sample  1 
Sample  2 


PMMA 


Sample  1 
Sample  2 
Sample  3 


Sample  1 
Semple  2 


PMMA 


Poly  Vinyl  Chloride 
( Ch loroethylene) 
Semple  1 
Sample  2 
Sample  3 
Sample  4 


Thermoelectric  Currente  from  Synthe.lzed  Polymer. 


Additive 


Son. 


Benzoyl  Peroxide 
(  .41  mole  %) 


( 1  mole  $) 


Benzene 

( 1  mole  %) 


Chi orobenzene 
( 1  mole  i) 


2-Methylpentene-l 
( 1  mole 


V  commercial; 


Semple 

Thickneas 

(cm) 


1,0 

0.* 


1.0 

0.4 


0.4 

0.4 

0.4 


0.4 

0.4 


1.0 
1.0 
0.4 
n  L 


Peak  Current 
( picoampa) 


O.kO 


3.6 

fc.3 


5.6 

4.6 


^.5 


O.53 

0-33 


32.5 

33-6 

56.0 
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4.  DISCUSSION 

The  disable  herein  i,  somewhat  speculative  since  th,  daU  u  ^ 

pert, .11),  complete.  Our  present  ideas  ere  briefly  strmarized  but  are 
subject  to  modification  based  on  future  results. 

a.  Model  of  Thermoelectric  Current  Generation 

The  results  on  the  synthesized  polymers  in  Table  lll-VI  are  quite 

revealing.  Pure  PMMA  with  no  additive  exhibits  a  very  small  current.  Few 

charges  appear  to  be  trapped  in  it  during  irradiation  at  35°C.  Thus 

either  or  both  positive  and  negative  charges  are  mobile  and  cause  recombina¬ 
tion. 

If  cci4  or  benzoyl  peroxide  are  added  to  PWA  the  observed  current 
increases  by  10  -  20  times  over  that  with  no  additive.  CC14  is  .,11  know, 

*°  react  with  electrons  by  dissociative  electron  attachment  i„  solid  organic 
matrioes.4  It  seems  likely  that  benzoyl  peroxide  can  act  in  a  simiLr  owner 
because  of  it,  weak  0-0  bonds  and  the  electron  affinity  of  oxygen. 
Alternatively,  it  is  probebl.  that  peroxide  and  oxide  radicals  from  benzoyl 
peroxide  are  produced  during  polymerizetion  and  can  act  as  electron  traps. 

It  should  further  be  pointed  out  that  electron  trap,  produced  by  dis.ocl.tiv. 
attachment  or  by  oxygen  confining  radicle  ere  quite  thermally  stable  and 

«ey  well  be  stable  above  reom  temperature  in  a  polymer  matrix.  That  i.  "d^" 
traps  are  produced. 

If  deep  electron  traps  are  provided  charge  is  stored  and  the  thermally 
generated  current  increases.  This  suggests  that  at  room  temperature  in 
pure  PMMA,  radiation-produced  electrons  are  mobile  and  recombine  with 
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trapped  holes  during  irradiation.  If  the  holes  were  also  mobile  in  the 
preaence  of  an  electron  trap,  recombination  could  still  occur  which  is 
contrary  to  our  observations.  This  tentative  conclusion  is  supported  by  the 
addition  of  2-methylpentene-1  (2M>-1)  to  PMMA.  In  solid  organic  matrices 
2HP-1  is  an  excellent  hole  trap  and  a  relatively  poor  electron  trap.4  PMMA 
doped  with  2MP-1  only  shows  a  twofold  increase  in  the  thermally  generated 
current.  A  much  larger  increase  (for  example  10-fold)  would  be  expected 
if  the  holes  were  mobile  and  the  electrons  were  trapped  in  pure  PMMA. 

The  twofold  increase  if  verified  in  additional  samples  is  more  likely  due 
to  the  weak  electron  trapping  nature  of  2MP-1. 

Both  benzene  and  chlorobenzene  act  as  poor  electron  and  hole  traps  in 
solid  organic  matrices.4  That  they  cause  no  increased  current  in  PMMA 
is  consistent  with  the  above  results. 

Polyvinyl  chloride  shows  a  very  large  current.  It  is  a  polyalkyl 
chloride  and  can  itself  react  with  electrons  by  dissociative  electron 

attachment.  Alkyl  chlorides  readily  show  this  reaction  in  the  solid 

1° 

phase.  In  contrast,  teflon  shows  negligible  current.  The  C-F  bond  is 
strong  enough  to  preclude  dis  ociate  electron  attachment.  (This  is  verified 
by  the  poor  response  of  linear  saturated  fluorocarbons  to  an  electron 
capture  detector  in  a  gas  chromatograph.)  Consequently,  the  results  on  PVC 
and  on  teflon  are  consistent  with  those  for  doped  PMMA. 

Commercial  PMMA  (Lucite)  is  thermally  polymerized  with  benzoyl  peroxide* 
thus  it  should  be  similar  to  PMMA  doped  with  benzoyl  peroxide.  This  has 
been  verified*  Lucite  shows  currents  of  similar  magnitude  to  those  of  PMMA 
doped  with  benzoyl  peroxide.  Thus  Lucite  is  convenient  to  work  with  since 
it  has  a  built-in  electron  trap.  It  should  b*  remembered  though  that  the 
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currants  observed  in  Lucite  are  characteristic  of  added  electron  traps  and 
not  of  pure  PMMA  itself. 

The  observation  of  thermally  generated  currents  indicates  that  one  or 
both  charge  carriers  are  mobilized  by  thermal  stimulation.  It  should  be 
noted  that  the  peak  current  occurs  near  65°C  which  is  well  below  the  glass 
transition  temperature  of  PMMA  at  90°C.5  Thus  the  temperature  of  detrapping 
seems  to  be  characteristic  of  the  trapping  site  rather  than  of  the 
unirradiated  polymer.  It  might  be  expected  that  if  the  electrons  were 
thermally  detrapped  that  the  temperature  at  which  the  peak  current  occurs 
should  depend  on  the  type  of  electron  trap.  However,  the  temperature  of  the 
peak  current  is  the  same  in  CCl+  and  benzoyl  peroxide  doped  PMMA.  This 
result  suggests  that  trapped  holes  are  mobilized  by  the  thermal  excitation. 
Hole  mobility  is  also  suggested  by  the  polarity  of  the  colder  electrode  as 
discussed  below. 

It  should  be  stressed  that  the  currents  studied  thus  far  are  thermo¬ 
electric;  they  require  a  temperature  gradient  for  their  generation.  Therefore, 
the  polarity  of  the  current  can  give  information  about  the  sign  of  the 
majority  charge  carrier.  The  existence  of  a  temperature  gradient  means  that 
there  is  a  net  propagation  of  lattice  vibration  waves  or  phonons  from  the 
hot  to  the  cold  electrode.  The  phonons  interact  with  the  mobile  charge 
carriers  and  cause  the  most  mobile  charge  carrier  to  drift  in  the  direction  of 
the  phonons.  Therefore  the  sign  of  the  cold  electrode  corresponds  to  the 
majority  or  most  mobile  charge  carrier.11'12  |n  all  measurements  on  samples 
thicker  than  0.10  cm  the  cold  electrode  has  been  positive  with  respect  to 
ground.  This  indicates  that  holes  are  the  mobile  species  and  supports  the 
suggestion  of  hole  mobility  based  on  the  additive  results. 
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b.  Stability  of  Trapped  Charge 

The  decay  curves  in  figure  1 1 1-16  can  be  divided  into  two  sections. 

The  initial  curved  section  is  formally  character jzed  by  a  reaction  order 
greater  than  one  while  the  straight  line  section  corresponds  to  first  order 
decay.  Similarly  shaped  decay  curves  have  been  observed  for  radiation- 
produced  hydrogen  atoms  in  solids.13  The  initial  curved  part  of  the  decay 
curve  can  arise  from  a  superposition  of  several  first  order  decay  curves 
with  different  activation  energies.  This  suggests  that  a  small  number 
of  charges  are  trapped  in  a  distribution  of  shallow  potential  wells.  Such 
shallow  traps  may,  in  fact,  account  for  the  small  currents  observed  in 
pure  PMMA.  The  majority  of  the  charges  are  trapped  in  deeper  potential 
wells.  These  sites  appear  to  be  characterized  by  a  single  trap  depth  or 
a  relatively  narrow  distribution  of  trap  depths  because  a  good  Arrhenius 
plot  is  obtained.  First  order  decay  corresponds  to  release  from  the  traps 
rather  than  to  charge  recombination  between  positive  and  negative  centers. 

The  decay  constants  derived  from  the  linear  sections  of  the  decay  curves 
show  a  good  fit  to  an  Arrhenius  plot  (Figure  111-17).  The  apparent 
activation  energy  is  20.7  kcal/mole  (0.9  ev)  in  Lucite  which  is  equivalent 
to  PMMA  doped  with  benzoyl  peroxide.  According  to  the  model  of  thermo¬ 
electric  currents  in  Lucite  postulated  above  this  activation  energy  is 
associated  with  the  hole  traps.  Experiments  are  in  progress  on  doped  PMMA 
to  determine  whether  the  activation  energy  is  independent  of  the  type  of  added 
electron  trap. 

The  mere  observation  of  a  large  activation  energy  indicates  that 
tunneling  is  not  the  mobility  mechanism.  At  present  a  thermally  activated 
hopping  mechanism  seems  likely  for  hole  mobility,  it  is  of  interest  to 
obtain  more  direct  information  on  the  mobility  mechanism. 


Sample  Thickness  Effects 
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c. 


Two  factors  could  affect  the  change  in  thermoelectric  current  with 
decreasing  sample  thickness,  the  temperature  gradient  and  the  trapped  charge 
distribution.  The  temperature  gradient  increases  with  a  decrease  in 
sample  thickness.  This  is  apparent  in  a  comparison  of  temperature  gradient, 
in  Lucite  samples  1.0  cm  thick  and  0.10  cm  thick  (Figures  MI-4  and  1 1 1-5). 
This  may  contribute  to  the  increase  in  current  when  the  sample  thickness  is 
decreased  from  1.0  cm  to  0.4  cm.  However  it  does  not  explain  the  change  in 
the  polarity  of  the  current  observed  during  measurement  of  samples  oetween 
0.10  cm  and  0.02  cm  thick  (Figure  111-14). 

The  charge  distribution  will  be  fairly  uniform  throughout  the  interior 
volume  of  the  thicker  samples.  At  the  surfaces  of  the  sample  the  distribu¬ 
tion  will  not  be  uniform  because  of  the  build-up  of  secondary  electron 
density  as  the  photon  flux  passes  f-om  a  lower  to  a  higher  atomic  number 
medium  (i.e.,  air  to  polymer).  !n  thicker  samples  the  effect  of  the 
nonuniform  charge  distribution  at  the  surfaces  is  small  conpared  to  the 
charge  distribution  throughout  the  volume.  As  the  sample  thickness  is 
decreased  the  surface  effects  are  accentuated.  More  work  is  necessary  before 
the  relative  importance  of  these  contributing  factors  can  be  ascertained. 

d.  Dose  Saturation  of  Traps 

Figure  II 1-9  indicates  that  the  normalized  total  charge  saturates 
after  a  dose  of  1.2  Mrad.  This  presumably  is  the  dose  at  which  all 
available  traps  are  filled  or  at  which  back  reactions  are  important.  It 
has  been  postulated  here  that  added  impurities  provide  the  trapping  sitee. 
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Anothar  possibility  is  that  free  radicals  formed  in  the  polymerization 
process  and  created  by  main  chain  scission  during  irradiation  could 
participate  in  charge  trapping.  Uhnishi  and  NittaU  have  measured  the  free 
radical  concentrations  in  PMMA  after  Co-60  irradiation.  They  report  that 
the  radical  concentration  saturates  at  a  dose  of  7  Mrad  and  is  a  linearly 
increasing  function  of  dose  to  about  2.5  Mrad.  Since  dose  saturation  of 
trapped  charge  occurs  at  much  lower  doses,  this  is  indirect  evidence  that 
free  radicals  do  not  provide  the  predominant  trapping  sites  in  PMMA.  The 
absence  of  large  currents  in  pure  PMMA  also  substantiates  this  conclusion. 
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IS  ABSTRACT 

Energetic  chemical  species  have  been  investigated  by  both  theoretical 
and  experimental  methods.  The  nature  of  the  chemical  bonding  in  Hj+  has 
been  studied  by  means  of  a  natural  spin  orbital  analysis  of  a  previously 
calculated  wavefunct ion.  Hj+  is  shown  to  rssemble  closaly  its  united  atom 
analog#  Li  ,  and  the  Hartree-Fock  energy  0f  Hj*  is  estimated  to  be  -1.301 
hartreos.  Also#  electron  density  plots  are  presented  in  order  to  clarify 
the  physical  picture  of  the  bonding  in  Ho*,  and  suggestions  are  made  aa  to 
how  further  calculations  might  proceed  most  officiently.  In  the  geminal 
theory  of  chemical  bonding  a  new  method  has  been  developed  for  solving  the 
coupled  integro-differential  equations  for  determination  of  tha  optical 
gsminals  for  a  system.  Garni nal  theory  has  also  been  generalized  to  include 
odd  electron  systems.  A  microwave  spectrometer  employing  Zeeman  modulation 


odd  electron  systems.  A  microwave  spectrometer  employing  Zeeman  modulation 
and  a  flow  call  has  bean  built  for  studying  transient  paramagnetic  species. 

,  Charge  storage  in  y-irradiated  polymethylmethacrylate  (PMAA)  at  room 
tsmpereture  has  been  studied  by  measuring  thermoelectric  currents.  By 
using  pure  PMUA  as  well  as  PMMA  doped  with  hole  and  electron  traps  it  has 
been  shown  that  added  electron  traps  but  not  hole  traps  lead  to  charge 
stabilization.  The  thermoelectric  currents  appear  to  be  due  to  hols 
mobility.  The  trapped  charge  in  doped  PMMA  decays  with  an  activation 
•nargy  of  0.9  sv.  Polyvinyl  chloride  stabilizes  even  larger  amounts  of 
charge  than  does  doped  PINA. 
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